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TECHNICAL NOTE D- 


OPTICAL PROBING OF SUPERSONIC 
AERODYNAMIC TURBULENCE WITH STATISTICAL 
CORRELATION. PHASE I: FEASIBILITY 

SUMMARY 


The theory and qualitative experimental results assessing the feasi- 
bility of measuring the statistical properties of supersonic turbulence by 
statistically correlating signals retrieved remotely with a laser schlieren 
system are presented. Also, some data retrieved with a laser shadow- 
correlation system are shown. Cross-correlograms and auto-correlograms 
computed on-line show that both of these systems can be used to retrieve 
flow-related signals sufficient for computing accurate and reproducible 
"peaks" of correlation. Additive tracers were not introduced because the 
schlieren and shadow-sensing modes were used. 

A statistical method for obtaining "one-shot" measurements of the 
decay history of turbulent structures in a stationary frame of reference is 
introduced, and results of practical applications of two types of these 
techniques are shown. The one-shot methods represent the only means by 
which the turbulence decay history can be computed from the same statistical 
sample of data (i.e. , the information can be computed from a single compo- 
site signal retrieved during one run of the facility). Theoretical analyses 
show that these one-shot techniques will also yield results in three- 
dimensional turbulent flow regimes with only minor modification of the beam 
arrangement (and in some cases with the insertion of time delays between 
signals) provided signal-to-noise ratios of the raw data time histories are 
not prohibitive. 

The application of parallel-beam geometry was used to increase the 
signal-to-noise ratio over that of crossed-beam geometry since only 
qualitative results were sought. Although the crossed-beam geometry should 
be used for quantitative measurements, the potential use of parallel beams for 
retrieving quantitative results from a restricted class of flows should be 
investigated. 


Electronically induced time delays are used as a means for (1) zoning 
the one-shot auto- and cross-correlograms, (2) avoiding peak overlapping, 
and (3) identifying the peak. 

With respect to the objectives of this investigation, the qualitative 
results presented in this report clearly indicate that optical remote probing 
of supersonic aerodynamic turbulent flows is feasible, without the use of 
tracers (using statistical correlation) . Further, correlation peaks computed 
from signals retrieved with either parallel- or crossed-beams during this 
investigation were flow-related, very reproducible, and readily identifiable. 

No discussion of spectra is presented and no attempt is made to 
experimentally establish what flow properties are measured. These should 
be the objects of a systematic test directed toward obtaining quantitative 
results. 


I. INTRODUCTION 


The majority of fluid flows encountered by the aerodynamicist are 
turbulent. Because of the complex nature of turbulence, the approach in 
the subsonic regime has been through the application of statistical methods 
to data obtained by solid probes inserted into the fluid. However, in super- 
sonic and hypersonic flows, these probes adversely influence the structure 
of the turbulence. One way to circumvent this difficulty would be to develop 
a reliable remote sensing tool for measuring the statistical properties of 
turbulence which would neither affect the flow field nor be adversely affected 
by it. The laser schlieren and laser shadow-correlation systems may, 
perhaps, represent the first generation of such a tool. 

The specific purpose of this publication is two-fold: (1) to present 

the theory and qualitative experimental results verifying the feasibility of 
optical remote sensing in supersonic turbulent flows employing the statistical 
correlation technique, and (2) to introduce and document the initial concept, 
theory, and reduction to practice of a statistical method which permits a 
one-shot measurement of the decay history of turbulent structures from one 
composite^ time history of flow information. 


1. The word "composite" is used in this report to imply a time history of 
random data composed of the algebraic sum of two or more random time 
histories of data. 
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The remote sensing tool used for this investigation employs two laser 
beams of light {X = 6328 A) which are retrieved by photodetectors after 
being influenced by the turbulent field. The time histories of these ac-coupled 
(fluctuating about a zero mean) signals are amplified and filtered, and 
statistical correlation methods are used to retrieve the desired flow informa- 
tion (e.g. , speed profiles, eddy lifetimes, turbulent length scales, and 
spectra) . 

The statistical one-shot autocorrelation technique constitutes methods 
whereby multiple time histories of random data are combined to form a 
single composite time history and an autocorrelogram is computed from the 
composite in such a way that a maximum amount of statistical flow informa- 
tion is retrieved with a minimum amount of time, equipment, and cost. 

The two basic types of one-shot correlations are referred to herein as the 
one-shot autocorrelation and the one-shot cross-correlation. Theoretical 
discussions, experimental results, and practical ramifications of these 
concepts are delineated in the main body of the report. 

The experimental objective of this work was to obtain accurate, 
reproducible, and readily identifiable correlation peaks with signals remotely 
retrieved from a supersonic turbulent flow that could be related to the most 
probable transit time of the disturbances. Previous attempts to obtain 
correlations in the supersonic regime have been hampered by (1) unknown 
anomalies in the data acquisition system (see Appendix C), and (2) the 
influence of facility-induced noise upon the optical system. 

The influence of facility- induced noise was reduced by utilization of 
the 7-inch Bisonic Wind Tunnel at MSFC, which operates at a very low noise 
level. This facility is described in detail in Appendix A. 

In addition to reducing the facility-induced noise, the power signal- 
to-noise ratio of the raw data was increased approximately an order of 
magnitude by placing the laser beams parallel in the two-dimensional 
turbulent boundary layer on a thin-plate model. The laser beams were 
parallel to one another in the horizontal plane and normal to the flow 
direction, thereby increasing the correlated signals between the beams. 

This correlation technique is analyzed theoretically in Section III, Paragraph 
A. 


Strong evidence that local information can be successfully retrieved 
from supersonic turbulent flows by using crossed beams is described in 
Section III, Paragraph C.4. Also, the cross-beam method was used to 


3 



investigate the effects the window boundary layers and interaction zones had 
upon the correlations computed from signals retrieved with parallel beams. 

It was found that, for the particular model design and thus the flow field, 
the boundary layers on the test section windows had no significant influence 
upon the correlograms. However, the interaction of the window and model 
boundary layers apparently dominated the measurements when the schlieren 
sensing mode was used in combination with parallel-beam geometry. Never- 
theless, this does not affect the conclusions of this investigation because the 
interaction zones near the windows are supersonic and turbulent. Also, the 
use of the laser shadow-correlation system reduced the contribution by the 
interaction zones approximately an order of magnitude. By crossing the 
beams in the vertical plane, these contributions to the correlograms can be 
avoided, as given in Section in. Paragraphs D. 4 and D. 5. 

The correlograms computed during this feasibility investigation 
represent our first encouraging measurements made in a supersonic turbu- 
lent flow with or without tracers and conclude the first phase of the wind 
tunnel cross-beam program. This has been conducted as an MSFC in-house 
research program with existing equipment and support. 


II. BACKGROUND 


In November 1968, a test was initiated to determine the feasibility 
of retrieving signals by optical remote sensing of supersonic turbulent flows 
for obtaining accurate and reproducible flow-related statistical correlations. 
The major problem in attempting to retrieve signals by this remote sensing 
technique has been the low signal-to-noise ratio. One of the major contri- 
butors to noise has been the facility- induced (mechanical and acoustical) 
excitation of the optical system and even of the flow field itself. Since the 
larger, more advanced wind tunnel and air- jet facilities being used in 
previous tests produced high noise levels, identifiable flow-related correla- 
tions could not be obtained in the supersonic regime. Because of the 
relatively small amount of noise produced by the Bisonic Wind Tunnel (BWT) 
facility at Marshall, this facility was selected for this investigation to alleviate 
the noise problem. 

The feasibility test in the BWT was planned in two parts, the objective 
of part one being to isolate the optical system from facility-induced noise. 

The possibility for successful isolation looked very promising during the 
initial investigations. It is difficult to accurately estimate the signal-to-noise 
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ratios obtained in previous facilities since facility-induced noise levels were 
always larger than flow-related signals, and also because it could not be 
determined that the correlations of signals were flow-related. 

In December 1968, a ratio of flow signal to noise of approximately 
10 was achieved in the BWT. A single laser beam was passed through the 
region of interaction of the shock wave and the free shear layer at position 1 
perpendicular to the flow (Fig. 1). The resulting signal is shown in Figure 
2A. The laser beam was moved to position 2 in the recirculatory region of 
the base. The signal at position 2 is shown in Figure 2B. All settings on 
instrumentation were the same for both runs. The experiment was repeated 
several times with the same results, thus implying that the increase in 
"signal" was flow-related^. The test objective of part one was achieved. 

The test objective of part two of the feasibility test was to obtain 
reproducible and readily identifiable correlation peaks related to the propaga- 
tion of disturbances in the two-dimensional supersonic turbulent boundary 
layer on a thin plate. An attempt to cross-correlate the signals from two 
parallel laser beams in the BWT was made during a demonstration on 
January 22, 1969. The laser beams were separated by 4. 5 inches and 
approximately 1/8 inch above the surface of the plate (Fig. 3). The signals 
from the two beams were amplified and cross-correlated in an analog 
correlator. The resulting cross-correlation, as a function of time delay, 
was displayed on a scope (Fig. 4A). The maximum correlation occurred at 
approximately 260 microseconds, which corresponds to a propagation speed 
of about 1440 fps^. The free-stream speed was approximately 1660 fps. 
Because of the nature of the run, a great deal of care in locating the beams 
in the flow was not taken. Therefore, the rheasurement revealed only that 
the peak occurred approximately at the expected time delay. 

Later investigations clearly showed that this correlation, though 
relatively weak compared to those computed later, was related to the 


2. Position 1 in Figure 1 was suspected to be one of considerable activity, 
because shadowgraphs of such shock-wave shear-layer interaction 
regions indicate this. 

3. The peak of the correlogram will be used to determine the speed of the 
disturbances. This is not necessarily correct because of the dependence 
of the correlation function on both space and time. In this report, 
qualitative results are sought, and therefore this simplification should 
be justified. 
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Figure 2A. On-line signal from single laser beam 
passed through the flow shown in Figure 1 — 
Beam at position 1. 
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Figure 2B. On-line signal from single laser beam 
passed through the flow shown in Figure 1 — 

Beam at position 2. 
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Figure 4A. Cross -correlogram from turbulent 
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Figure 4B. Cross-correlogram in free shear layer 
of 2-D wedge (M = 2.0). 


propagation of turbulence in the boundary layer on the plate and in the inter- 
action zone of the plate and window boundary layers. Further, a considerable 
amount of water vapor was present that extinguished the laser beams by 
scattering (schlieren effects were also present). Nevertheless, this very 
crude measurement represented the first encouraging result from our attempts 
to retrieve signals from supersonic turbulent flows. 

In November 1968, during part one of the feasibility test, an attempt 
had been made to cross-correlate the signals from one laser beam which had 
been split with a beam splitter into two beams of equal power (a second laser 
was not available at that time). The purpose of this run was to determine if 
there were correlated harmonic signals present in the flow like those obtained 
during the preceding test in the 14-inch Trisonic Wind Tunnel (see Appendix 
B). The two beams were passes through the free shear layer of the wedge 
shown in Figure 1. The correlation at zero time lag (Fig. 4B) was due 
partially to the correlated laser noise. The large peak at about 100 micro- 
seconds gives a speed of 1080 fps for the beam separation of 1. 3 inches. 

The available supporting equipment and the test objective of part one would 
not allow the pursuit of this interesting result. Also, the facility-induced 
noise had not been fully investigated. Thus, no significant level of confidence 
could be put in the probability that the correlation or any portion of it was 
actually flow-related. In retrospect, it seems reasonable to suspect that 
there was a significant flow-related contribution to the correlogram. The 
indirect influence of this particular run upon the test results is, perhaps, of 
interest and will be discussed later. 

After the cross-correlation was obtained in phase two (Fig. 4A), it 
was necessary to determine the origin of the signals from which it was 
computed. It was evident from the shape of the subsequent cross- 
correlograms of the individual beams (Figs. 5A and 5B) that there was little 
or no periodicity present in the signals. Furthermore, because of the near- 
zero value of these cross-correlations at zero time delaj and the absence of 
correlated periodic noise, it could be deduced that the correlated facility- 
induced noise was very small. Finally, the shape of the autocorrelograms 
indicated that the correlated signals were of the wideband variety [ 1 ] . Two 
additional conditions had to be met to establish the fact that the correlations 
were truly flow-related: (1) The peak had to occur at a time delay which 
corresponded to the approximate expected speed of the disturbances. 
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and (2) condition (1) had to hold as the beam separation (|) was varied. 
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Figure 5A. Autocorrelation of signal 1. 
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Figure 5B. Autocorrelation of signal 2. 
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The first runs to check these necessary conditions were made with the 
upstream beam (beam 1) and the downstream beam (beam 2) located in the 
wake of the model at positions 3 and 4, respectively (Fig. 3). The resulting 
cross-correlogram of the signals is shown in Figure 6A. The beam separa- 
tion and the time delay corresponding to maximum positive correlation gave a 
speed of 1415 fps, which was at least reasonably close to what was to be 
expected (free-stream speed was approximately 1660 fps). Then, beam 2 
was moved upstream such that the beam separation (|) was decreased by 1 
inch. The cross-correlogram for this case is shown in Figure 6B. The 
maximum correlation corresponded to a speed of 1409 fps. Also, as would 
be expected, the strength of the peak correlation increased. The only 
difference between these runs, other than beam separation, was that the 
amplitude scale (voltage) on the oscilloscope was changed for the second run 
to accommodate the increase in amplitude of the correlation peak. 

These results, verified by later experimental data discussed in the 
following sections, provide reasonable experimental evidence supporting the 
acceptance of the flow-related nature of the correlation peaks. 

Several new experimental and theoretical techniques were employed 
during this investigation which are: (1) application of paraUel beams for 

remote probing of a two-dimensional turbulent flow, (2) the laser schlieren 
and laser shadow-correlation modes of retrieving signals from a supersonic 
turbulent flow, (3) the one-shot auto- and cross-correlation methods for 
statistical analyses of random time histories, and (4) the theory of induced 
time delays for peak identification purposes. 

Perhaps the most significant product of this feasibility test should be 
attributed to the attempt to correlate the random fluctuations of two beams of 
light originating from a common source. In itself, the splitting of a laser 
beam is, of course, common practice, but not previously applied using the 
cross-beam method. Nevertheless, the use of parallel beams from a common 
source led to the initial one-shot autocorrelation concept described in Section 
III, Paragraph B. 1. , and then to the one-shot cross-correlation, peak 
identification by the method of electronically-induced time delay and a visual 
on-line display of the one-shot auto- and cross-correlograms. These methods 
are discussed in later sections. 


12 



Cross - Correlation Cross - Correlation 

(Volts)^ _ (Volts)^ 



Time Delay (Milliseconds) 



Time Delay (Milliseconds) 


Figure 6B. Cross -correlogram from measurement 
made in wake of plate with 1 inch less beam separation 
than Figure 6A, 
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IN. DISCUSSION 


A. Remote Sensing by Optical Correlation 
with Parallel Laser Beams 

A theory of remote sensing of clear air turbulence using parallel 
laser beams (or other parallel light sources) in the supersonic regime is 
discussed, and experimental results are presented. 


1. THEORETICAL DESCRIPTION OF A LASER SCHLIEREN SYSTEM 

In Figure 7, the beams from two lasers, Lj and L 2 , are directed 
through the two-dimensional turbulent boundary layer shown in Figure 3. 

The flow is supersonic (M = 2. 0), and the boundary layer is fully devel- 
oped along the thin plate. The beams, which are parallel to one another and 
perpendicular to the flow direction, are separated in the flow direction 
(along the x-axis) in such a manner that the eddies which intersect beam 
LjDj at time, t, also intersect beam L 2 D 2 at a later time, t +t . 

Since the flow is approximately two-dimensional, the time-averaged 
flow properties should be approximately equal across the test section (i. e. , 
as measured along either laser beam). This does not necessarily require 
that the instantaneous flow properties across the test section be equal. 

The turbulence passing through one of these beams exhibits random 
fluctuating localized gradients of flow properties; e. g. , local density gradient. 
As the eddies pass through a particular location on the beam, the resulting 
fluctuations of flow properties produce fluctuations of the local Poynting 
vector^, s, of the laser beam. 

If ]3(y,t) represents the local instantaneous angular deviation, as 
measured from a horizontal reference line perpendicular to the flow direc- 
tion, then 

s • f = ! si cos I3{y,t) , (1) 


4. The Poynting vector is a vector tangent to the path of the beam with 
magnitude equal to the beam power. 
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where f is a unit vector along the y-axis. Rearranging equation (1) results 
in 


/3(y,t) =cos"^ 



( 2 ) 


This fluctuation, /3(y,t), is caused by a fluctuation of the local index of 
refraction, n(x,y,z,t), and is related to the local fluctuation in the density 
gradient normal to the Poynting vector, s, for a particular wavelength of 
light. The local index of refraction, n(x,y,z,t), along the beam is pro- 
portional to the local density, p(x,y,z,t). Provided the index of refraction 
is close to unity, as it is for most gases, the empirical Gladstone-Dale 
equation [3] is applicable. This equation, an empirical relationship between 
the index of refraction and the density of gas, is of the form 


^ — — = constant = a (3) 

Pi P2 


so that 


n^(y,t) =1+0! Pj^(y,t) (4) 

for any point along the beam. 

Let us assume that n and p are represented by 

n(x,y,z,t) = n(x,y,z) + n'(x,y,z,t) (5) 


and 


p(x,y,z,t) = p(x,y,z) + p'(x,y,z,t) 


( 6 ) 
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where n' and p' are the fluctuations about the time-averaged values of 
n and p. Substituting equations (5) and (6) into (4) gives 

n=n+n'=l+Q!p+ ap' (7) 


and 

n' = cup' . (8) 

It is desired to relate the beam deflection angle, /3 (y,t) , to the 
component of the index of refraction gradient which is perpendicular to the 
Poynting vector (path of beam). Since the index of refraction for air is 
very close to unity and since i3(y,t) will be small, 

<»(y.t) = ig-dy , (9) 

Equation (9), the derivation of which can be found in most textbook pre- 
sentations on the schlieren method of flow visualization [2], represents 
the fundamental relationship upon which the schlieren method is based. 
Substituting equation (7) into (9) yields 

c*3(y,t) i 1^1 

Since n - 1 « 1 (n for .air is approximately 1. 00027), equation (10) can 
be simplified one step further: 

The deflection of the Poynting vector away from the y-axis at a particular 
location on the beam (17) will be obtained by integrating along the beam 
from y = -i to y =r}: 


9n 3n' 

dx dx 


dy 


( 10 ) 
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( 12 ) 


/3(i7,t) = / dy + / I?- dy 


-i 


ax 


The first integral on the right-hand side of equation (12) represents the 
temporal-average angular beam deflection at t) , and the second integral 
represents the fluctuation of the beam deflection about the temporal average. 
Thus, 


/3(T],t) = |8 (t]) + /3’(^?,t) 


(13) 


Substituting equation (7) into (12) provides a relationship between /3 and 
p: 


/3(r?,t) 


a 


/ 


-a 


ax 


dy + O' 


/ 

-i 


ax 


dy 


(14) 


The beam deflections in these relationships are in the x-y plane only. It 
will be shown later why the deflections in the z-direction can be neglected. 
Also, another assumption will be introduced: the magnitude of the time- 
averaged beam deflection, ^ (i ) , is such that the knifercdge (which is used 
to produce the fluctuations in current monitored by the photodiode) allows 
one-half of the laser beam light to pass into the photodiode when the beam 
deflection equals j8(i). Thus, equation (14) reduces to 


/ If dy 


(15) 


This assumption is convenient and can be applied experimentally. The 
fluctuations of /3' can be related to the fluctuating output signal, i(t), of 
a photodiode, as shown in Fignre 7. 

A view taken along the laser beam (Fig. 8) shows the knife-edge, 
the eye of the photodiode, and the laser beam cross section. The beam is 
shown in a deflected position and in its time-averaged position (centered). 
Let 


I 


d 


time averaged current output of the detector power supply 
when knife-edge is removed. 
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Figure 8. Effect of knife-edge on the deflected beam. 

A = cross-sectional area of laser beam, 

D = diameter of laser beam. 

A = beam displacement perpendicular to and measured from 
the knife-edge. 

AA = the change in area of the beam cut by a knife-edge due to 
a deflection of the beam off center. 

Then, Figure 8 shows that 

AA = ^ sin-1 (2A/D) + ^ n/ - 4A^ . (16) 

For small A , 

AA = DA . (17) 
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The photodetector output per unit area of the laser beam is 



A 7tD^ 


( 18 ) 


To obtain the fluctuation of the output voltage (or current) from the 
photodetector, i, we now multiply the output per unit area by the decrease 
(or increase) in beam area cut by the knife-edge due to a beam displacement 
of A: 


41 A 

i(t) = — , for (A « D) 

TTD 

Further, A is related to )3'(i,t): 


(19) 


A = d • /?' (i,t) 


( 20 ) 


Substituting equation (20) into (19) gives 

■ji3'(i,t) (21) 

The group of terms in parentheses on the right-hand side of equation (21) 
is a measure of the sensitivity, s, of the system [4], 

4df 

- - -J- ■ 

The sensitivity of the photodiode with respect to beam position was 
assumed constant, as was the intensity profile across the laser beam 
diameter. These assumptions, however, may be far from true and should 
be checked. By placing a lens between the knife-edge and the photodiode, 
the beam movement on the photodiode can be considerably reduced. 

Substituting equation (22) into (21) gives 



i(x, t) = s • /3' (x, i, t) 


(23) 


20 





2 . 


CROSS-CORRELATION OF SIGNALS 


Equation (23) holds for both beams so that 


ij(x,t) = Si/3|(x, i,t) (24) 

i 2 (x + |,t) = S 2 • /32 (x + ^,i,t) , (25) 


Substituting equation (15) into equations (24) and (25) provides relationships 
between i and dp'/dx, 


- s.. / dy. 

-i ° 


i2(x + ^,t) = S2<a / ■^-^2S2L±i±IiAll 


(26) 


(27) 


which hold along a particular eddy '’streamline. ” 

The signal, i 2 , from the downstream beam, is delayed by r seconds 
and the time-average magnitude of the cross product, i^ • i 2 , is calculated. 
This represents the cross -correlation of the signals and I 2 [I = l/2f^ + i(t)], 

or the cross-covariance of the fluctuating components, i^ and i 2 . The 
cross-correlation of i^ and i 2 will be used in the following discussion. 

The cross-correlation, R(^,r ), is defined as [1] 

1 ^ 

R(^,t) = lim — f ii(x,yi,t) • i 2 (x + |,Y 2 , t + T ) dt , (28) 

T — « o 


where T is the averaging time. It is assumed that T is large enough so 
that 


R(^,t) = Y 


T 

/ ii(x,yi,t) 
o 


+ |,y2,t + t) dt 


(29) 
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is a good assumption. Also, stationarity is assumed (i. e. , time-averaged 
statistical values are independent of time). 


In equation (29) , R(|,t) will be maximum when 4 is delayed by 
the most-probable transit time of the common disturbances (t= t^), 

provided the signal-to-noise ratios and averaging time, T, are compatible 
and will allow the desired peak to rise well out of the correlated noise floor. 
Otherwise, the correlated noise will influence the shape and maximum 
position (t^) of the peak. Correlated noise can often hide the peak 

altogether. The latter is likely to occur when the peak height is approxi- 
mately equal to, or less than, the level of the correlated noise floor. 


Substituting equations (26) and (27) into (29) gives 


1 r 

B(i.r) = ~ / 


!,£< / 


8Pi'(x,yi,t) 


ax 


dyi 


s^Q, J l£2iE±liM±ll 


dx 


dy2 


dt 


(30) 


Rearranging the order of integration, we obtain 




where the bar indicates the time-average of the product inside the braces. 
Thus, equation (31) provides one interpretation of what the cross-correlation 
represents. 

The crossrcorrelation of the signals ij and i, from two parallel 
beams is equal to a constant (siS 20 i^) multiplied by the double line integral, 
one taken along each beam, of the two-point-product time-average value of 
the fluctuating density gradient component in the flow direction. 

Another, and perhaps a better, description is the one which shows 
the sequence of physical relationships between the parameters being measured 
experimentally. 
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Beginning with the on-line analog correlator and moving opposite to 
the flow of data, the following physical events are taking place (see Appendix 
A, Fig. A-11); 

(a) The scope displays the cross-correlation function versus time 
delay (t ) which is received from the analog correlator. 

(b) The correlator is computing the cross-correlation from the two 
electrical signals, ij and i 2 , originating from the photodiodes. 

(c) The time history of the signals are directly proportional to the 
angular deflections of the laser beams, respectively: 


i(t) = s/3' (i, t) 

(d) The angular deflection /3' (i ,t) of each beam is proportional 
to the fluctuating component of the density gradient in the direction of flow, 
where the mean is an instantaneous spatial average taken along the laser 
beam; i. e. , from equation (15): 


/?’ U,t) 


a 


i 


/ 


9p'(x,y.t) 

dx 


dy 


(15) 


Multiplying and dividing the right-hand side of equation (15) by the width of 
the test section (2i = L) gives 


/3’ (i ,t) =2Ha 


f 9p'(x,y,t) 

2i K dx ^ 


Thus, at any particular instant (tp), 

fap'(x,y,t)' 

£l 


/3 (f ,t ) =2f O' 
P 


ax 


avg. over L 


(23) 


(32) 


The substitution of t = t is made only to emphasize that the average is 

P 

taken along the beam as if the flow were frozen in time. This is not implying 
that the beam-length average value of dp'/dx does not vary with time. 
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because it does. This results from the speed of light being much greater 
than the flow speed. Thus, the effect is (or can be) of considerable 
advantage to two-dimensional aerodynamic remote-sensing with statistical 
correlation. The advantage is a decrease in the integration time required 
to retrieve correlated statistical information from data in the presence of 
noise because of the automatic optical integration performed by the laser 
beam. For the two-dimensional case, the signals ii(t) and i 2 (t) represent 
the fluctuating component of the density gradient, (dp'/dx), spatially 
averaged along the respective beams. Therefore, the correlator is operating 
on data that are representative of instantaneous averages of random samples 
of data. The convergence, therefore, will be faster than for the three- 
dimensional case where the beams must be crossed. Further, and much 
more important, the signal-to-noise ratio is at least an order of magnitude 
greater for two-dimensional flow fields similar in flow rate and geometrical 
size. This is not meant to imply that optical remote sensing of three- 
dimensional flows is not possible. 

From the experimental point of view, these comments can be 
summarized as follows: 

(a) The correlation is representative of the similarity in the two 
electronic signals, ij and i 2 , received by the correlator after they have 
been influenced by the filters, amplifiers, and all other physical connections 
between the photo power supply output and the correlator input (see Appendix 
A, Fig. A- 11). If the correlator is performing properly, this is simply a 
matter of experimental fact. 

(b) The signals are proportional to the angular deflections of the 
respective beams about a time-average value. 

(c) The signal represents the fluctuating (beam-average) x-component 
of the density gradient. 


B. The One-Shot Autocorrelation 

1. THE CONCEPT 

In Figure 9, a single laser beam is passed normal to the flow and 
along the particular station of interest. After passing through the test 
section, the beam is reflected downstream from mirror Ml to mirror M2. 
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Figure 9. The one-shot autocorrelation for determination 
of convection speed in two-dimensional turbulence. 

From mirror M2, it is reflected back across the flow at a predetermined 
location and parallel to the first pass through the flow. A photodetector, 

Dl, receives the beam after its two passes through the test section. The 
signal is amplified and the analog correlator computes the autocorrelation 
versus time delay (autocorrelogram) , as shown at the bottom of Figure 9. 

The autocorrelogram has two predominant peaks, one at zero time 
delay, as expected, and the second at a time delay, t^, corresponding to 

the most-probable transit time of disturbances between the two positions 
where the beam passes through the test section. Although the beam arrange- 
ment is similar to that of the parallel beam case, the desired flow informa- 
tion is contained in a single signal rather than in two separate ones. The 
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second peak occurs at on the autocorrelogram because the signal 

received by the photodetector, Dl, contains the same (or similar) informa- 
tion twice and at different times. The most-probable separation in time 
corresponds to t 

m 

This statistical technique can be used in analyzing random data such 
as those produced in turbulent fluid flows. The technique, referred to in 
this report as the "one-shot" autocorrelation, is described as follows: 

• The autocorrelogram of a single composite time history 
composed of the sum of two or more statistically cor- 
related random signals which sufficiently lag one 
another will exhibit a correlation peak for each possible 
pair of the signals. 

2. MULTIPLE SIGNALS 

The extension of the concept would be to reflect the beam through the 
flow three or, perhaps, four times, as shown in Figure 10. Although this 
technique will work in practice, it is preferable to avoid such a sensitive 
optical arrangement. 

In Figures llA and IIB, an equivalent arrangement is shown. One 
laser beam is split into four beams of equal intensity and directed through 
the flow. Each beam has a separate detector that receives the signal. The 
signals are added and the resulting signal is equivalent to the original single 
beam case. It is not necessary to use a single light source. Each detector 
could have an independent light source, or two or more detectors may share 
light from the same source. 

The one-shot autocorrelation technique can, in theory, be applied 
to two- or three-dimensional flows. The only change required is a 
rearrangement of the beams so that no two are parallel and so that the 
effective correlated volume that is common to the beams is small enough to 
be acceptable. However, a longer integration time will be required for 
three-dimensional flows, as compared to two-dimensional flows, because of 
the lower signal-to-noise ratio of the former. 


3. PEAK IDENTIFICATION 

Consider a single-time history of bror.d-band random data, F(t) 
(i.e. , distributed over a wide range of frequencies) which is composed of 
two random-time histories fj(t) and f 2 (t) such that 
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Figure 10. Multiple pass parallel beam arrangement for the one-shot autocorrelation. 
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Figure llA. Schematic of multiple beam arrangement 
with separate detectors. 

One-Shot“ Autocorrelation 
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F(t) =fj(t) +f2(t) 


(33) 


and the mean values of the time histories are zero. The autocorrelation of 
F(t) is 


A(t) = lim f [fjCt) + f 2 (t)] [fj(t + T ) + f 2 (t + T )] dt (34) 
1 — 2 T — o 


A(t) = fj(t) • fj(t + T) +fj(t) • f2(t + r) +f2(t) • fj(t + T) 
1—2 


+ f2 (t) • f2 (t + t) (35) 


A(t) = Aj(t) +Ri2(t) +R2i(t) +A2(t) . (36) 

1 — 2 


Inequation (36), A and R represent autocorrelation and cross- 
correlation, respectively. The subscripts indicate the signals involved in 
the particular correlation, as well as the order of the cross-correlation 
operation with regard to the time histories. 

Since, 


R2i(t) = Rj2(-t) 


(37) 


Substituting equation (37) into (36) gives 


A(t) = Aj(t) +Rj2(t) +Rj2(-t) +A2(t) . (38) 

1—2 


If we assume, for the present, that fj(t) and f 2 (t) are statistically 
independent and thus not correlated, equation (38) reduces to 


( 39 ) 


A(t ) = Aj(t) + A2(t) 
1—2 


29 



For this case, the one-shot autocorrelogram will have only one large peak® 
at T = 0 similar to the autocorrelogram shown in Figure 12A. The 
magnitude of the peak at t = 0 is 


A(0) = Ai(0) +A2(0) . (40) 

1 — 2 


Beginning again with equation (38), we assume that fj(t) and 
f 2 (t) are identically equal for all t when f 2 (t) is delayed by seconds 

and compared to fi(t); i. e. , at 


T = T 


m 


fi(t) = f 2 (t + T ), for all t 

i. y 


(41) 

(42) 


where t is the time lag of maximum correlation. For this case, 
m 


F(t) = fi(t) + f2(t) 


(43) 


and substituting equation (42) into equation (35) gives 


A(t ) = Aj(t ) +Ri2(t ) +Rj2(-t ) +A2(t ) 
m ^ m m m m 

1—2 


(44) 


If we assume that is positive and large enough to prevent an interference 

between the peak at and the zero time lag peak on the autocorrelogram, 

5. Peak implies the maximum positive peak in a local region on the auto- 
correlogram. The natural shape of the auto- and cross-correlograms, 
near the time lag, , generally exhibits three or more peaks depend- 
ing upon the frequency range of the correlated portion of the signals. 


30 


Autocorrelation of Signals Autocorrelation of Signol Two 

Qne and Two (Volts)^ 1 , (Volts)^ 


I 



Time Delay (Milliseconds) 

12A. Autocorrelation of signal 2 versus time delay. 



Time Delay (Milliseconds) 

Figure 12B. Autocorrelation of signals 1 and 2 
versus time delay. 
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we can deduce from equation (44) that 


since 


A(t ) = Ri2(t ) 
m m 

1 — 2 

(45) 

Aj(t ) = 0 
^ m 

(45a) 

A2(t ) = 0 
^ m 

(45b) 

R<2(-t ) = 0 

(45c) 


Equation (45) states that, if fj(t) andf 2 (t) are statistically random-time 

histories and have maximum correlation when f 2 (t) is delayed in time by 

T seconds, then the cross-correlation of two signals is equal to the one- 
m 

shot autocorrelation provided that is sufficiently large. Therefore, 

there should be no problem in identifying the peak at because there 

should be only one (Fig. 12B). There is, however, a need to look further 
into the composition of the one-shot autocorrelogram. 

Consider the case of the addition of three random signals where 


F(t) = fj(t) +f2(t) +fg(t) 

The one-shot autocorrelation of F(t) is 


(46) 


A(t) = Ai(t) +A2(t) +A3(t) +Rj2(t) +Rj3(t) +R23(t) +R2i(t) 
1—3 


+ R32(t) +R3i(t) . (47) 
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Using the general form of equation (37), we obtain 


A(t) = A^(t) +A 2 (t) +A 3 (t) +Ri 2 ('’’) +R23("r) 

1—3 


+ Ri 2 (-t) +R23(-t) +Rig(-T) 


( 48 ) 


If it is assumed that all peaks resulting from the cross-correlations in 
equation (48) lie in the positive time-lag range only, the last three cross- 
correlations in equation (48) drop out. 


A(t) = Ai(t) +A 2 (t) +A 3 (t) +Rj 2 (t) +Ri 3 (t) +R23(t) . (49) 

1—3 


Again, let us assume that the time lags between regions of correlation on 
the one-shot autocorrelogram are large enough to prevent overlapping of 
peaks. Let Tjg be equal to the time lag of maximum correlation between 
fj(t) and fglt) , ti 3 the time lag of peak correlation between fj(t) and 
fg(t) , and T 23 the time lag of maximum correlation between f 2 (t) and 
fg(t). Then, the value of the one-shot autocorrelation at the above 

will be 


A(ti2) 

1—3 

= Ri2('^12) +1^23 ('^ 12 ) 

(50) 

A ( Ti3 ) 

~ ^13(T13) I^23^'^13) 

(51) 

1 — 3 



A ( T23 ) 

= R23^'’’23) +Ri2("^23) ^^13 ^ '’’23 ) 

(52) 


1 — 3 


If we assume that no overlapping of regions of correlation occur, then 
equations (50), (51), and (52) reduce to 
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A ( T 12 ) - 12 ( 12 ) 

1—3 


A(ti3) = Ri3(ti3) (54) 

1 — 3 


A(t23) - 1^23^ '>'23^ (55) 

1—3 


Figure 13 shows a three-beam arrangement, where the parallel beams pass 
over the plate and through the turbulent boundary layer. The free-stream 
flow speed was approximately 1660 fps. The signals from the two photo- 
detectors, A and B, were added, amplified, and filtered. The resulting 
signal was autocorrelated on-line. The autocorrelogram is shown in the 
top photo of Figure 14. The peak at the zero time lag is characteristic of 
all autocorrelograms of the wideband time history. The peak at r = 0. 128 ms 
is apparently a result of the correlation between beams 1 and 2. It will be 
shown later that this peak is actually a result of the correlation between beams 
1 and 2 summed with that between beams 1 and 3, The be.am separations, 

(i. e. , Ii2 distance between beams (1) and (2), etc. ) were 


4 i 2 =2.63 inches and ^^3 =4.95 inches 
which give average speeds of 

and 


Ui3 


1 13 
12ti3 


1525 fps 
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Figure 13 . Parallel beam arrangement with three laser beams. 

(Both optical and electrical addition of signals are used. ) 

The beams were approximately 0 . 2 inch above the plate, near the outer 
region of the boundary layer. The speed Uj2 hardly seems reasonable 
because the free-stream speed was 1660 fps; the speed Uj3 is closer to the 
expected value. However, there should have been a third peak representing 
the correlation between the second and third beams, 1^23(^23). Dividing the 
beam separation, I23 . the speed, U13 , should represent a good approxi- 
mation to the location of the peak (T23) (^23 = ^13 “ ^12 = 2.32 inches). 
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Figure 14. One-shot autocorrelations as obtained with a three- 
parallel beam arrangement as shown in Figure 13. 
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T 23 = _ — = 0. 127 ms 

12Ui3 

It appears that the peak at 0. 128 ms is actually a result of the overlapping of 
two peaks, and thus does not represent the independent location of either. 

To check the above logic, the first beam was moved 1 inch downstream. This 
reduced the beam separation distances, ^13 > 1.63 and 3. 95 inches 

respectively. However, the separation between the second and third beams, 

^23 , was not changed. The "one-shot" autocorrelogram for this beam 
geometry is shown in the bottom photo of Figure 14. 

From Figure 12B, 


Tj 2 = 0 . 082 ms 
T 13 = 0. 203 ms 
T 23 = 0. 123 ms 


The corresponding speeds are 

Uj 2 = 1660 fps 
Ui 3 = 1521 fps 
U 23 = 1570 fps 


The purpose of this run was to show that the peak at T 23 did exist. 
The value of U 12 should not be considered accurate because the peak at 
T 12 = 0 . 123 ms is too close to T 12 to be considered accurate. There are, 
however, at least two ways by which peak overlapping can be avoided. These 
will be discussed in Section III, Paragraph B. 5. It is important that the 
number of peaks expected to appear for a particular number of beams is 
known in advance, as well as their expected location on the one-shot auto- 
correlogram. 

In general, the number of peaks on the one-shcft autocorrelogram 

will be 
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No. of peaks 


1 + 


(56) 


m-1 

Yj (m - k) 
k=l 


provided that none are destroyed by overlapping. The peak at zero time lag 
has been included in equation (56). 

The one-shot autocorrelation could be used to study the character- 
istics of the decay of turbulent structures from data obtained from a single 
run of a wind tunnel or jet facility. The peaks resulting from correlations 
between the first beam and each of the downstream beams give a picture of 
the lifetime of the average eddy in a stationary frame of reference. In 
addition, the peaks resulting from correlations between the second beam 
and the downstream beams represent the eddy lifetime a second time, but 
beginning at a different location (i.e. , ^2 further downstream). Likewise, 
the series of peak associated with the third beam and each beam downstream 
from it provides still another history of the decay. Therefore, the one-shot 
autocorrelogram is composed of at least one decay history (two beams) and, 
where more than two beams are used, there will be m-1 of such histories, 
each history beginning a distance downstream from the last equal to the beam 
separation (Fig. 15). 

If the assumption is made that the decay history is independent of 
position over the distance between the first and last beams, we see that all 
of these time histories fall on the curve of the first time history, which 
begins with the first beam (Fig. 16). It follows, for this case, that a peak 
resulting from correlation between two downstream beams provides the same 
information as the correlation between the first beam and an additional beam 
having the same beam separation as the two downstream beams. 

This characteristic of the one-shot autocorrelation is important with 
respect to the number of beams required to determine the decay history. 

For example, three beams provide four points on the decay curve, and four 
beams provide seven points. 


4. INTERPRETATION OF THE ZERO- TIME LAG PEAK 

The zero time-lag value of the decay curve is not represented by the 
zero time-lag value of the one-shot autocorrelation. However, this value 
can be calculated to a good approximation in many cases. 
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Figure 15, The one-shot autocorrelation of a three-time 
history composite showing decay envelopes of 
turbulent structure. 

Consider the general case of m number of time histories where 
f Jt) represents the ith time history and 

m 

F(t) = Yj . (57) 

i=l ^ 


The general form of the one-shot autocorrelation is 

T m m 

A(t) = Urn / Z f-(t) • E f.(t + T)dt . (58) 

1 — m T — «> o i=l j=l ^ 
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Equation (58) can be represented more conveniently in terms of auto- and 
cross-correlations, as was done previously. 


m m 

A(t) = Yj Tj • <59) 

l^m i=l j=l •' 


For the applications considered here, 




for X > Y , and 

Rxx<^) = A^(t) (61) 


where X and Y are particular values of i and j . Therefore, 


m mm 

A(t) = Y A (t) + Y Y Rii<T) , j>i . (62) 

l-*m 1=1 1=1 j=2 •' 


Equation (62) is the general form of the one-shot autocorrelation for this 
particular type of application. For the present, assume that the beams have 
equal sensitivities, s, and that the zero time-lag value of the individual 
autocorrelations of the f.(t)^s are equal. Then, 

m 

A(0) = Y A.(0) = mA^(O) (63) 


and 

Ax<0) = ^ [A(0)1 . (64) 

l-*m 


|| 
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If it is further assumed that the signal-to-noise ratio approaches infinity as 

the beam separation approaches zero between the first beam and any one of 

the downstream beams (i. e. , no noise when =0), then the value of the 

1 ] 

one-shot autocorrelation is equal to the cross-correlation at t = 0. 


Ax(0) ^XY 


(65) 


For this particular case, the curve representing the trace of peaks as a 
function of t on the one-shot autocorrelogram (i.e. , the curve representing 
the rate of turbulence decay) will pass through the point 


T = 0, 


- A(0) 
m . 

1 — m 


This value, at r = 0, is an upper limit and must decrease as the signal- 
to-noise ratio decreases, because, at t = 0 on any autocorrelogram, the 
noise is correlated as strong as the signal. 


5. TIME DELAYS AND ZONING OF THE ONE-SHOT ' 

AUTOCOERELOGRAM 

Overlapping of regions of correlation on the one-shot autocorrelogram 
prohibits certain combinations of beam separations and limits the minimum 
separation between any two of the beams. There is a solution to the over- 
lapping problem referred to herein as the method of induced time delay. 
Electronic time-delay techniques are employed in many fields. In this 
report, there are three purposes of the induced time delay: (1) to avoid 

peak interference on the one-shot correlogram, (2) to identify a correlation 
peak partially lost in correlated noise, and (3) to zone the one-shot auto- or 
cross-correlogram, thereby providing means for identifying particular pairs 
of signals associated with a particular correlation peak on a one-shot auto- 
or cross-correlogram. 

With respect to (1) above, Figure 17A shows a one-shot auto- 
correlogram where the measurement was made with two parallel laser 
beams separated by 1. 90 inches in Mach 2. 0 flow. The beams were located 
in the turbulent boundary layer. 
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2 One - Shot 

Cross Correlation (Volts) § Autocorrelation (Volts) 



Time Delay (Milliseconds) 


17A. One-shot autocorrelation of two signals added. 



Time Delay (Milliseconds) 

Figure 17B. Cross-correlation of signal 1 with signal 2. 
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The effect of overlapping is clearly evident by comparing Figures 17A 
and 17B. In Figure 17B, a cross-correlation was made between the two 
beams which eliminate the overlapping problem when only two beams were 
used. 


There is a solution to the overlapping problem. If the second signal 
of Figure 17A had been delayed by electronically introducing a known time 
delay, At, the peak at ~ ^06 ms would have occurred at 

T ' = T + At. If a time delay of 0. 100 ms had been introduced into the 
mm 

signal from the downstream beam, the one-shot autocorrelogram would 
have looked like the one shown in Figure 18, which is a much more accurate 
correlogram from the standpoint of studying the nonzero time-delay correla- 
tion. 


In particular, induced time delays provide flexibility to the one-shot 
autocorrelation technique and extend its range of application. For example, 
the speed and dire . ‘ n of atmospheric ground winds can be theoretically 
determined by the one shot autocorrelation method. 

Figure 19 is a schematic diagram of an atmospheric wind detection 
system composed of six signals retrieved by optical remote sensing 
techniques®. Five of the six detectors are arranges in a fan to determine 
which signal from the fan is correlated strongest with the signal from an 
”up-wind" location. It is also desired to determine the time lag of maximum 
correlation. Ideally, this can be accomplished with the one-shot autocorrela- 
tion method using induced time delays of known value. 

As shown in Figure 19, the signals from the fan B, C, D, E, and F 
are each delayed electronically and then added together. The resulting 
signal is added to signal A. The one-shot autocorrelogram of the summed 
signal should, theoretically, appear as shown at the bottom of Figure 19. 

Induced time delays can be used to zone the one-shot autocorrelogram 
for peak identification. If no time delays are introduced, the peaks would 
probably overlap. Even if no overlapping occurs, there would be doubt as 
to which signal from the fan produced a particular peak. 


6. Memo by B. H. Funk, Jr. , entitled "Theoretical Application of the 'One- 
Shot' Autocorrelation Technique to the Determination of Atmospheric 
Ground Wind Speed and Direction," Memo No. R-AERO-A-68-37, 

Dated December 11, 1968. 
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The induced time delay, At, is large compared to the expected 
transit time, Since signals B, C, D, E, and F are delayed iDy At, 

2At, 3At, 4At, and 5At seconds, respectively, a peak resulting from 
correlation between signals A and B would be expected to occur at 
T '= T + At: likewise, for signals A and C, t ' = t h- 2At, etc. Also, 

the induced time-delay method provides a graphical relationship between 
wind direction and the abscissa of the one-shot autocorrelogram. It follows 
that by comparing the one-shot autocorrelograms of the same data but with 
different induced time delays, peak identification should be enhanced. 

The induced time-delay method of peak identification should apply 
just as well to cross-correlation, perhaps better. 


6. NOISE ANALYSIS^ 

The major disadvantage of the one-shot autocorrelation is that the 
noise-to-signal ratio, taken with respect to a particular time history con- 
tained in the composite, is greater than the noise-to-signal ratio of the 
particular time history taken separately. This is necessarily true, because, 
in a region on the autocorrelogram where two signals of a composite time 
history are correlated, all other signal and noise components represent 
noise in the computation. For applications where the noise-to-signal ratios 
of the individual time histories are large, the one-shot integration time may 
become very long to achieve the same accuracy as that obtained from the 
individual cross-correlations. However, this has not been proven experi- 
mentally. 

In the following, let 

n = noise-to-signal ratio of composite time history, 
a 

n = noise-to-signal ratio of individual time history. 

m = number of time histories composing the composite 
time history. 

i = instantaneous value of the composite time history, 
c 

7. The conventional signal-to-noise ratio has been inverted here because 
the mathematics are simpler. 
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instantaneous value of the kth individual time history 
(k = 1,2, . . . , m). 


instantaneous signal of i. . 

k 


i = instantaneous noise of i . 

n,, k 


a. Assumptions 

(1) The noise-to-signal ratio of the individual time histories are 


equal. 


(2) The rms value of the individual time histories are equal. 

(3) There is statistical stationarity. 

Derivation 

From the definitions, it follows that 


m 

i = ii + io + i-5 + . . . + i = ) , i. 

c 12 3 m . k 

k=l 


( 66 ) 


1“ = rms value of composite time history, 

c 


ic = ^ / (ii + i2 + • • • dt 

a o 


(67) 


i^ = ij + i| + . . 
c ^ ^ 


. + 2iji2 + 


( 68 ) 


The cross-product terms approach zero as T -► «=. Thus, 

a 


_ _ m 

i2 :2 . j2 , , ,2 _ V ;2 

'c - ^ '2 + . . • + 1,^ - L \ 

k=l 


(69) 
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If the noise-to-signal ratio is defined as the ratio of the rms noise in the 
composite to the rms value of the signal, then 


n = 
a 



1/2 


and 


n = 
a 




1/2 


( 70 ) 


(71) 


If the individual noise-to-signal ratio is defined as the ratio of the rms 
noise to, the rms signal, taken with respect to the individual time histories, 
then 


n = 



(72) 


and 


n2 = 



(73) 


Now substituting equation. (69) into (71), 


n 

a 



i! 





By assumption 2, 


(74) 



(75) 
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(76) 






= 7 


and 


= . . . = 

Sj nj S2 H2 s n 

By assumptions 1 and 2, 


(77) 



1 = r = . . . - i 

Hi n 2 n 

Substitution of equation (76) into (69) gives 


(78) 

(79) 


i^ = m i^ 
c 


(80) 


and from equation (77) , 




s 



Replacing i^ in equation (80) by the expression in (81) yields 


(81) 


i^ = m< 


(i^ + i2 \ 

\ s n/ 


(82) 


Substitute equation (82) into (71); 



( 83 ) 


50 



Substitute equation (73) into (83): 


n = 
a 


I m ( 1 + n^ ) 



( 84 ) 


Equation (84) is the rms noise-to-signal ratio of a composite time history 
taken with respect to one of the individual signals and expressed as a function 
of the rms noise-to-signal ratio of the individual signal. 

The power noise/signal ratio of the composite time history is 


n 

pa 


n‘ 


The power noise-to-signal ratio of an individual time history is 


(84a) 


n 

P 


n 


(84b) 


Then, from equation (84), 


n=m(l+n)-l . ( 84c) 

pa p 


In Figure 20, n 

3 , 

For large values of n, 


is plotted versus 


m for several values of n. 


n = \l m n 
a 


(85) 


For the range of values plotted in Figure 20, equation (85) is not a bad 
approximation even for n = 4, m = 2. For n = 0, which corresponds to 
the case where there is no noise in the individual time histories, 


n = \Tm - T 
a 


( 86 ) 
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Ratio of Noise - to - Signal of 
0 Composite Time History, n 


Number of Time Histories in Composite, m 


[ C ' 

rms Values And Noise ■ tO“ Signol Ratios of The 
Individual Time Histories Are Assumed Equal. 


Figure 20. Noise-to-signal ratio of a composite versus 
number of time histories added. 




It is seen from equation (86) that the minimum noise-to-signal ratio of the 
composite time history is unity. 

The practical limitation is not set by the increase in noise-to-signal 
ratio directly; rather, it is the increase in integration time required to 
obtain the same results that would be obtained by calculating the individual 
cross-correlations with equal confidence levels. The increase in integration 
time can be described as the ratio of T to T. Here, T is the integration 

time required to compute the one-shot autocorrelation and T is the integra- 
tion time required to compute a cross -correlation between two of the individ- 
ual time histories. The relationship is® 


T 2 . 

a _ 2 m^ - 1 

T ~ ^ ~ (n® + 1)® + 1 


(87) 


In Figure 21, the functional relationship between T /T, m, and n 
is shown. From this figure the following observations can be made; 

( 1 ) The minimum value of T is 

3 


T = (2.5)T 
a 


which corresponds to the case of adding two identical signals to form the 
composite, i. e. , n = 0 and m = 2. However, for this case, T is also a 
minimum. Therefore, adding two, three, or four such signals should be 
considered practical. Verification of this has already been presented for 
m = 2 and m = 3. 

(2) The percentage of increase in integration time due to adding 
time histories with noise-to-signal ratios n= 3 is essentially the same for 
fixed m; i. e. , to a good approximation. 


8. Memo by M. Y. Su, entitled "Statistical Error and Required Averaging 
Time for Evaluating Cross-Correlation of Cross-Beam Data," 
Nortronics -Huntsville Memo No. M-794-8-362, dated May 24, 1968. 
The n used in the memo is the square of the n used here. 
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T 

, for n = 3 
T 

(3) The limits T /T are 

H. 


+ 1 
2 



The increase in the required integration time, then, represents a very 
important limitation upon the practical application of the one-shot auto- 
correlation method. There is, however, another limitation that is not so 
critical. This is the increase in computation time. 

The ratio of computation times is given by 
T 

V = ( 

^a (number of peaks) • T 


As explained in Section III, Paragraph B. 3. , the number of peaks on the one- 
shot autocorrelogram is 

m 

number of peaks = 1 + ^ (m - k) . (56) 

k=l 


The ratio of computation times is 


^a = 


m 


1 + Yj (m - k) 
k=l 


Thus, 


r„ = 


m(m - 1) + 2 T 


(89) 


(90) 
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which has upper and lower limits; 


(m^ + 1) ^ < 2m^ 

m(m-l)+2 “ ^a ^ m(m-l)+2 


( 91 ) 


In Figure 22, y is plotted versus m. This figure represents the case 

9 . 

where all peaks of the one-shot autocorrelogram are considered as useful 
data. 


C. The One-Shot Cross-Correlation 

1. STATEMENT OF THE CONCEPT 

The second kind of one-shot correlation technique is the one-shot 
cross-correlation, which is described in the following statement: 

• The cross -correlation of a single time history 
with a composite time history, the single time 
history containing a random signal that is 
statistically correlated to two or more random 
signals contained in the composite, will result 
in a cross-correlogram exhibiting peaks of 
cross-correlation between the signal in the 
separate time history and each of the signals 
in the composite, the signals in the composite 
sufficiently lagging one another in time. 

Two important differences between the one-shot autocorrelation and 
the one-shot cross-correlation are as follows: (1) The one-shot cross- 
correlogram does not have the large peak at zero time delay representing the 
mean square value of the composite signal, and (2) the peaks resulting from 
correlations between all possible pairs, other than those pairs formed 
between the signal in the separate time history and each of those in the 
composite, are not present on the one-shot cross-correlogram. 

These differences represent a loss of information. There are cases, 
however, where it is desirable to eliminate this additional computation. 

For example, when the noise-to-signal ratio of the individual time histories 
of a composite used in computing a one-shot autocorrelation are large. 
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Figure 22. Ratio of computation time (7^) versus 
number of time histories in composite (m). 
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usually the error in evaluating the zero time-delay point on the eddy decay 
curve is large. Also, the required integration time may be very long. In 
such a case, a considerable advantage can be gained by using the one-shot 
cross-correlation. The advantages are (1) lower noise-to-signal ratio and, 
therefore, a shorter required integration time, and (2) the capability to 
compute the cross-correlation peak at zero time delay directly without an 
electrically induced time delay, which is necessary with the one-shot 
autocorrelation. 


2. MATHEMATICAL DEVELOPMENT 

In the following, consider the cross-correlation of a single time 
history, with a composite time history, F(t). F(t) is composed 

of ni individual time histories such that 


m 

F(t) = X . (92) 

k=l ^ 


Assume that each time history in F(t) and f^(t) represents random time 

histories and that each can be expressed as the sum of signal plus noise 
where a signal in F(t) represents that portion of a time history which is 
correlated to the signal of Then, 


f (t) = f (t) +f (t) 
o s n 

o o 


F(t) 


m 


Z 


k=l 


f (t) +f (t) 
s, n 

k k 


Furthermore, assume for each fj^(t) that 


f^U) 


= lim 

T ^ oo 


T 


T 

o 


(93) 


(94) 


(95) 
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and 


f (t) = lim — r f (t) dt = 0 
o T o 

T -» 00 o 


(9G) 


The one-shot cross-correlation of f (t) with F(t) is defined as 

o 


1 r 

R(t) = lim f f (t) • F(t + r) dt 

1 — m T 


T 

oo CO 


(97) 


and 


c m 

E(r) = lim f f (t) • E 

1— m c o k=l 

c 


(98) 


It is assumed that T is sufficiently long, such that 

c 


m 


R(t) = ^ f f^(t) • Z + dt 
1 — m CO k=l 

is a good approximation to equation (98). 

Substitute equations (93) and (94) into (99): 


(99) 


T 

R(t) = / 

f (t) +f (t) 

m 

• z 

f (t + t) + f (t + t) ] 

T 

1 — m CO 

CQ 

O 

0 

1 

k=l 

1 

ca 

1 


dt 


( 100 ) 

Performing the indicated multiplication and taking the summation outside 
the integral, we obtain 
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E(t) 

1 — m 


Z 


k=l 


_1 

T 


c 


T 

/ fg 

o o 


f (t + x) dt + 
®k 


T 

/° 

o o 


f (t + r) dt 
n. 


T 

c 

+ r f (t) 

^ n 
o o 


f (t + r) dt + 

\ 


T 

/ 

o o 


• f (t + t) 
n, 



( 101 ) 


The last three integrals on the right-hand side of equation (101) are 

approximately zero because the time-average value of the products of 

signals with noise, and noise with noise are zero; i. e, , the noise component 

of f (t) is not correlated with the signal and noise components in the 
o 

composite, F(t). 

Equation (101) reduces to 

T 

m c 

R(t) = Yj ^ S fg (t) • fg (t + t) dt . (102) 

1 -*• m k=l coo k 


The summation in equation (102) represents the sum of the cross-correlations 
of the single time history f^(f) with each of the time histories in the com- 
posite evaluated at the same time delay, t . 


R(t) 

1 — m 


Z 


k=l 




(103) 


where the subscripts, ( ) , , indicate the signals cross-correlated, as 

ok 

well as the order of the cross-correlation. 

If it is assumed that the signals in the composite are introduced 
sequentially such that each signal in time lags the previous signal sufficiently 
to avoid peak interference, the one-shot cross-correlation expressed in 
equation (103) becomes 
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(104) 


R(t 
1 — 



m 


R At .) 
oj oj 


where is the time delay corresponding to the peak of the cross- 

correlation R At .) between f (t) and f. (t + r) in the composite. All 

OJ OJ o j 

other cross-correlations in equation (103) evaluated at are zero by 

the assumption that overlapping does not occur. 


m 


Z 

k=l 


k^j 


R , (t .) = 0 
ok OJ 


(105) 


Therefore, the one-shot cross-correlogram will be composed of as many 
peaks as there are time histories in the composite, each of which contains 
a random signal statistically correlated to the random signal in f^(t)- 


number of peaks = m. 


(106) 


The time delay between the kth and (kth + 1) peaks is t , ^ - t , 

o,k+l o,k 

which is sufficient to avoid overlapping. 


3. EXPERIMENTAL VERIFICATION 

In Figure 23, three laser beams, parallel to each other and normal 
to the mean flow direction, are directed through the turbulent boundary layer 
on a thin plate. The free-stream Mach number is approximately 2. 0. The 
beams form a plane that differs only slightly from a plane parallel to the 
surface of the plate. Beam 1 from the upstream laser passes through the 
boundary layer, and a knife-edge blocks approximately 50 percent of the 
laser light, preventing it from reaching the photodiode. The uninterrupted 
portion of the beam is monitored by the photodetector A . The fluctuations of 
the laser beam perpendicular to the knife-edge causes fluctuations in the 
amount of light monitored by the photodiode that converts these fluctuations 
into an ac signal. 
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The downstream laser beam is split into two beams of equal power. 
Beam 3 passes directly through the boundary layer. Beam 2 is reflected 
upstream and then through the boundary layer at a location that is 2. 32 inches 
upstream of beam 3 and 2. 63 inches downstream from beam 1. After beam 2 
passes through the boundary layer, it is reflected downstream. Beams 2 and 
3 are added optically by monitoring both beams with the same photodetector 
B. The fluctuations of beams 2 and 3 are retrieved by use of a knife-edge in 
the manner described from beam 1. 

The cross -correlation of the time history monitored by detector A, 
with the composite time history monitored by detector B, is shown in Figure 
23 (lower right). The two peaks predicted by the one-shot cross-correlation 
theory are clearly evident and provide strong experimental support to its 
validity. 

To compare the one-shot cross-correlogram obtained above with the 
one-shot autocorrelogram calculated from a composite of the same three 
time histories, the output of detector A was electronically added to the output 
of detector B, and the autocorrelation of the composite was computed as a 
function of time delay. This one-shot autocorrelogram is shown in Figure 23 
(upper right). The structure of this one-shot autocorrelogram was previously 
discussed in Section III, Paragraph B.3. 

In Figure 24, a one-shot cross-correlation is shown with three beams 
in the composite time histories. During this particular run, the second beam 
in the composite was not positioned on the knife-edge. Instead, it was located 
on the edge of the photodiode opposite to the knife-edge (in the downstream 
direction). Thus, the edge of the photodiode performed the same function as 
the knife-edge, except the sign of the signal was reversed. This inverted the 
second peak on the one-shot cross-correlogram and demonstrates the 
schlieren principle very well. 


4. NOISE ANALYSIS 

For the same individual time histories, the integration time required 
to compute the one-shot cross- correlation would be less than the integration 
time required to compute the one-shot autocorrelation, both having the same 
number of peaks on their correlograms (excluding the peak at zero time 
delay on the one-shot autocorrelogram). 
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Figure 24. One-shot cross-correlogram with three time histories 
in the composite. (Sign of second signal was negative. ) 


In the following, 


T 

c 

n 


‘o<‘> 

i (t) 
c 

m 

i^(t) 


= integration time for a one-shot cross-correlation. 

= rms noise/signal ratio of a signal in a composite time 
history. 

= instantaneous value of the independent time history. 

= Instantaneous value of the composite time history. 

= number of beams in the composite time history. 

= instantaneous value of the kth individual time history 
in the composite. 


T = integration time required to compute the cross-correlation 

between the independent time history and one of the time 
histories in the composite. 
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The rnis noise-to-signal ratio of i is defined as 

o 



(107) 


and from Section III, Paragraph B, the noise-to-signal ratio for a signal in 
a composite is 


n = [m(l + n^) - 1] 
c 


(108) 


In equation (108), it has been assumed that the noise-to-signal ratio, n, 
and the mean-squared values of the individual time histories in the composite 
are equal. Also, statistical stationarity is assumed. 

The ratio of integration times is 


T n^ n^ + n^ + n^ + 2 

c c c 


T n*^ + 2n^ + 2 


(109) 


It has been assumed in equation (109) that the noise-to-signal ratio and 
mean-squared value of the independent time history, i^ , are equal to those 

of an individual time history in the composite; i. e. , 


n = n 
o 


o 



Rearranging equation (109) yields 


T 


c 


T 


(n^ + 1) 

(n^ + 1)^ + 1 



( 110 ) 

( 111 ) 


( 112 ) 
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Substituting equation (108) into (112) and rearranging gives 


T 


c 


T 


(n^ + 1) ^ 

(n^ + 1) + 1 


(m - 1) 


(113) 


or 


T 

c 


T 


(ni - 1) 
(n^ + 1)2 + 1 


(114) 


The upper limit of equation (114) is 


lim (T /T) -*m , (115) 

c 
n 

and the lower limit is 


lim (T /T) - , (116) 

n-*0 ^ ^ 


or 


m + 1 ^ 

2 



(117) 


Equation (114) shows that for n= 3 the ratio of integration times 

(T /T) is, to a very good approximation, equal to the number of individual 
c 

time histories in the composite: 

T 

= m; for n = 3 . (118) 

In Figure 25, T /T is plotted versus m for various values of n. 
c 
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Figure 25. Integration time ratio versus m for a 
one-shot cross-correlation. 


From equation (56), the number of peaks on a one-shot autocor- 
relogram are (excluding the peak at r = 0) 

m 


E - k) 

k=l 


The number of peaks on the one-shot cross -correlogram are 


P = m 
c c 


However, for comparison purposes, it should be required that 


number of peaks = 


a 


= ^ (m2 - m ) 
2 a a 


(119) 


( 120 ) 


( 121 ) 


P = P 
c a 


( 122 ) 


Substituting equations (120) and (121) into (122) gives 


m = - (m2 - m ) 
c 2 a a 


(123) 


Replacing m in equation (87) with m , we obtain 

a 


= m‘ - 


m - 1 
a 

(n2 + 1)2 + 1 


(124) 


which is the one-shot autocorrelation integration time ratio. Replacing m 

in equation (114) with m , we obtain 

c 


T m - 1 

c _ ^ c 

T “ ”^c “ (n2 + 1)2 + 1 


(125) 
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Dividing equation (124) by (125) yields 


T (n^ + 1)^ + 1 

a _ a 

T “ m (n^ + 1)2 + 1 
c c 

and substituting equation (123) into (126) yields 


(126) 


2m2 (n2 + 1)2 + 2 

3 . 

m (m - 1) (n2 + 1)2 + 2 
a a 


From equation (127), the following conclusions can be drawn: 

(1) The lower limit of the integration time ratio is 

2m2 + 2 
a 


(127) 


lim (T /T ) . - 

a c m(m-l)+2 

a a 

n — 0 

(2) The upper limit is 
2m ^ 

lim (T /T ) 
a c 


(128) 


n — 00 


m - 1 
a 


(129) 


(3) The integration time required for the one-shot autocorrelation is 

always longer than that required for the one-shot cross-correlation. The 

minimum value occurs in the limit as m and is equal to 2. 

a 


a 


c/m in 


lim (T /T ) — 2 
a c 

m — ► oo 
a 


(4) The maximum value of T /T corresponds to the limit as 

a 

n — 0 for m =3 and is equal to 3. 
a 
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lim (T /T ) — 3 
^ a c 
n -- 0 

m =3 
a 


(5) For comparison purposes, the minimum value of m is 3 

a 

because the one-shot autocorrelogram for m -2 has only one peak, other 

a 

than that at t = 0, and the corresponding one-shot cross-correlation reduces 

to the standard cross-correlation. By equation (123), m =1 for m =2. 

c a 

These equalities imply that there is only one time history in the composite of 
the one-shot cross-correlation. 


(6) The upper limit is approached rapidly as n increases. For 


n ^ 3, 


2m 


m - 1 
a 


n g 3 


(130) 


with less than 1 percent error. 


(7) The range of T /T is 
a c 


T 

2. 50 ^ ^ 3. 00 

c 


for all combinations of n's and m 's between the limits defined below. 

a 


0 ^ n ^ °o 


3 ^ m ^ °o 
a 

Therefore, the required integration time for a one-shot autocorrelation is 
250 to 300 percent longer than that required for a one-shot cross-correlation. 
This statement is, of course, restricted to the assumptions made above. 
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In Figure 26, T /T is plotted versus m for several values of n. 
a c a 

If the power noise-to-signal ratio had been used in the analysis in place of the 

rms noise-to-signal ratio, conclusion (3) would have been the only one 

affected (n ^ 3 would have been n ^ 9 since n = n^). 

P P 

The ratio of the computation time required for the one-shot cross- 
correlation, , to that required for m number of separate cross- 
correlations (one for each peak on the one-shot correlogram) is 


y 


c 



m - 1 

m[ (n^ + 1)^ + 1] 


(131) 


The lower limit is 


, . / V m -H 1 , , 

lim (y ) — (132) 

c 2m 

n — 0 

and the upper limit is 


lim(y^)— 1 , (133) 

n — oo 


so that 


^1, m^2 (134) 

2m ^c 

Thus, from equation (134), the computation time required to compute the 
one-shot cross-correlation is always less than the time required to corhpute 
the m number of cross-correlations separately. However, as n increases, 
the ratio, y ^ , rapidly approaches its upper limit. For n ^ 3, 


y^ = 1.0, for n = 3 


(135) 
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Number of Beams In Composite of 
"One -Shot" Autocorrelation, m^ 


Figure 26. One-shot auto- and cross-correlation 
integration time comparison. 

with less than 1 percent error. Therefore, the one-shot computation time 
is never a disadvantage. Actually, does not accomit for computer 

printout time and other data-handling procedures minimized by the one-shot 
method. 

Figure 27 shows plotted versus m for n = 0, 1, and °° , 
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D. Cross-Beam Discussion 


1. GENERAL DISCUSSION 

The "cross-beam" theory [5] is essentially applicable to the 
measurements discussed below, which were made using crossed beams. 
However, as described in Section III, Paragraph A. 1. , the signals were 
retrieved from the flow by employing the schlieren principle. The cross- 
beam runs differ from the parallel beam runs in that the correlated signals 
from crossed beams are retrieved from a "localized" region defined by the 
segments of the two beams through which the same disturbances pass. The 
length of the beam segment is approximately the size of the average dis- 
turbance passing through and common to both beams. Because the beams 
are crossed, the flow disturbances that do not pass through both beams, 
and thus do not cause correlated fluctuations in the signals, represent 
noise. As a result, the cross-beam method should require more integration 
time than the parallel beam method. Nevertheless, the crossed-beam 
geometry offers many more advantages than disadvantages. The development 
of the cross-beam method is the ultimate objective. 

In the following discussion, a crossed-beam geometry was used to 
investigate the contribution to the correlations computed using parallel beams 
resulting from (1) the boundary layer on the test section windows, (2) the 
interaction of the wall boundary layer with the boundary layer on the plate, 

(3) the wake from the indentations in the plate running parallel to the 
windows (Fig. A-8), and (4) a combination of these. Also, a crossed-beam 
geometry was used to make measurements in the wake of the model. 


2. EFFECT OF THE BOUNDARY LAYER ON THE TEST SECTION 

WINDOWS 

In Section III, Paragraph A. 1. , it was assumed that the statistical 
properties of the turbulence along each of the laser beams were constant. 
This assumption is not valid for positions near the test-section windows 
because of the influence of the wall boundary layer. However, the magnitude 
of the influence was the important point to be determined. 

The necessity for this investigation results from the use of parallel 
beams. Because the average transit time of the disturbances between two 
parallel beams decreases from the window to the outer edge of the boundary 
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layer, the effect should produce an unsymmetrical contribution to the cor- 

relogram distributed from a time delay (t) slightly larger than the transit 

time corresponding to free-stream speed between the beams to larger values 

of T. Also, the total contribution from the wall boundary layer should be 

small compared with the correlation resulting from the disturbances along 

the beams that are greater than the boundary layer thickness, 6 , away 

w 

from the windows (Fig. 28). 

Two runs were made using parallel beams separated by 1. 87 inches. 
The first run was made with the beams passing through the boundary layer 
on the thin-plate model (Fig. 3). The centerlines of the beams were 0. 05 
inch above the surface of the plate. The cross -correlation of the two 
signals is shown in Figure 29A. The time delay of maximum correlation, 

, is equal to 0. 112 ms (0. 106 ms plus 0. 006 ms due to parallax). The 

same result was obtained when the run was repeated. 

The beams were elevated to approximately 1. 07 inches above the 
plate. One run was made and repeated twice. From the shadowgraph of the 
flow (Fig. 3), it can be seen that this is well above the major portion of the 
boundary layer (6 « 0. 3 inch). 

All settings on the instrumentation were the same for all three runs. 
The correlogram of one of the runs is shown in Figure 29B. The corre- 
lograms of the other two runs were essentially the same as the first. 

A comparison of Figures 29A and 29B would lead to the conclusion 
that the boundary layer on the test- section windows through which the parallel 
beams pass does not significantly affect the correlogram. The observation 
that the amplitude of the signals decreased by approximately a factor of 5 
when the beams were in the elevated position well above the boundary layer 
supports this conclusion. 

Although the conclusion is reasonable, we cannot assume that the 
wall boundary layer 1. 07 inches above the plate has the same effect on the 
correlogram as the flow resulting from the interaction of the wall and plate 
boundary layers near the surface of the plate. This region of interaction 
may contribute heavily to the correlation; therefore, it was necessary to 
design a test that would separate the correlations along the beams, such that 
the contributions caused by the interaction zone could be separated from those 
caused by the boundary layer on the model. 
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Figure 29A. Cross -correlation of separated parallel beams 
0. 05 inch above plate. 
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Figure 29B. Cross -correlation of separated parallel beams, 
1. 03 inch above plate. 


3 . 


INVESTIGATION OF THE BOUNDARY LAYER INTERACTION 
ZONE 


Figure 30 shows the plan view of the test section. The two laser 
beams were crossed such that the point of intersection was in the middle of 
the test section and the plane formed by the beams was approximately parallel 
to the surface of the plate. 

Two runs were made with this geometry and the cross -correlograms 
were computed on-line. Because of the geometry of the beams and the flow, 
the cross-correlogram is symmetrical about the origin (t = 0). The 
positive time-lag range was computed during the first run (Fig. 30A), and 
the negative range was computed during the s econd run (Fig. 30B). Cross- 
correlations between the beam segments AC and BC correspond to the 
positive range of the cross-correlogram, and cross-correlations between 
beam segments DC and EC correspond to the negative range. 

Two dominant peaks occurred on the cross-correlogram, one at 
+ 0. 167 ms and the other at - 0. 167 ms ( + 0. 158 ms +0. 009 ms due to parallax) . 
The common disturbances producing the correlations represented by these 
large peaks can be estimated by multiplying the approximate speed of transit, 
measured with parallel beams at the same distance above the plate (1510 fps) 
by the transit time (time delay) taken from the cross-correlogram. 

I = U • T 

approx. m 

= (1510) (12) (0. 167) 10~2 

= 3. 02 inches. 

Thus, the beam segments through which the common disturbances causing 
these large peaks pass are separated by approximately 3. 0 inches. From 
the beam geometry, the beam segments are shown very near the windows 
in the interaction zone of the wall and plate boundary layers. 

To check this result, the position of beam intersection was changed 
(Fig. 31A). For this run, the beams were crossed in the wall boundary 
layer on the detector side of the tunnel. The plane formed by the beams 
was essentially the same as that of the previous case ( 0. 20 inch above the 
surface of the plate). 

The cross-correlogram of beam 1 with beam 2 representing the cross- 
correlation between beam segnaents AC and CB is shown in Figure 31B. The 
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Figure 30A. Cross-correlation measurements 
of the influence of the interaction zone upon parallel 
beam correlations — Positive time-lag range. 


Figure SOB. Cross-correlation measurements 
of the influence of the interaction zone upon parallel 
beam correlations — Negative time-lag range. 












time lag of maximum correlation is 0.263 ms (0. 248 + 0. 015 ms for parallax). 
Again, using the mean speed of disturbances that was measured with parallel 
beams, the corresponding distance of separation can be estimated. 


^ =U-T = (1510)(0.263) 10-3(12) 
0S u* m 

= 4. 77 inches. 


Thus, the beam segments through which the common disturbances producing 
the correlation peak pass are separated by approximately 4. 8 inches. From 
the geometry of Figure 31A, the beam segments shown are separated by 
4. 8 inches, as measured in the direction of flow, and are well inside the 
interaction zone near the window. 

The positive correlation at zero time lag in Figure 3 IB was expected 
and is caused by the correlation of the disturbances passing through the 
intersection of the beams. 

Figures 32A and 32B show the negative time-lag range of the cross- 
correlogram representing cross-correlation of disturbances passing through 
beam segments DB and EB. The zero time-lag peak was expected; however, 
the large peak previously obtained in the negative time-lag range with the 
beam intersection in the center of the test section was not computed. An 
acceptable explanation for this is not presently available. 

For the next two runs, the horizontal plane of the beams was 
elevated 1. 55 inches above the surface of the plate. No other changes were 
made for these two runs that were made to evaluate the contribution of the 
wall boundary layer to the correlogram. 

Figure 33A shows the positive time-lag range, and Figure 33B shows 
the negative range. These cross-correlograms corroborate the results 
obtained using parallel beams; i. e. , that the wall boundary layer alone would 
have very little influence upon the correlograms computed from signals 
retrieved from the plate boundary layer if the two flows did not interact. 

Finally, it can be concluded, but not without an element of uncertainty, 
that the interaction zone near the window and plate model is a region of intense 
turbulence and contributes significantly to the magnitude and shape of the 
correlograms computed from signals retrieved with parallel laser beams 
presented in this report. 
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Figure 32A. Schematic of beam geometry — Negative time-lag range. 
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Figure 32B. Negative time-lag range of cross-correlogram 
for beam geometry shown in Figure 32A. 
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33A. Positive time-lag range of cross-correlogram 
for beam geometry shown in Figure 32A. 
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Figure 33B. Negative time-lag range of cross-correlogram 
for beam geometry shown in Figure 33A, 
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In general, these results do not invalidate measurements taken during 
this test series, since the test objectives were primarily to ascertain the 
feasibility of applying electro-optical remote sensing techniques to making 
measurements of supersonic turbulent flows, and to demonstrate the one- 
shot techniques without extensive alterations to the tunnel or expenditures of 
time and manpower. From this standpoint, the primary test objectives were 
achieved, since parallel beam correlation represents an optically integrated 
ensemble average correlation of eddy transit times across the entire test 
section. However, before quantitative measurements are made using parallel 
beams, a two-dimensional flow field with minimized edge effects should be 
sought, possibly through improvement in model design and sensing mode. 

Other attempts to minimize edge effects were also examined. One 
method investigated is a type of shadow-correlation technique that is 
described in Section III, Paragraph D. 5. Crossing the beams in a vertical 
plane that is perpendicular to the flow also effectively eliminates edge 
effects, but the cross-beam method is rather cumbersome for measuring 
the velocity profile or other flow properties near flat plates because of 
geometrical limitations. Transparent models would, perhaps, allow cross- 
beam measurements, but passing a beam through a transparent model 
introduces a new problem of beam stability that should be investigated. 
Nevertheless, cross-beam measurements should be quite successful and 
convenient for making measurements in the wake. The next paragraph 
discusses one such attempt to obtain these cross-beam measurements in the 
turbulent wake of the thin-plate model under investigation. 


4. CROSS-BEAM MEASUREMENT IN TURBULENT WAKE 

This paragraph discusses the feasibility of making cross-beam 
measurements in the turbulent wake of the thin-plate model. More space is 
available in the wake for geometric beam arrangements than in the boundary 
layer on the model. However, the beams could not be crossed perpendicular 
to one another [5] in a vertical plane perpendicular to the flow because the 
BWT does not have windows available for passing a vertical beam through 
the flow. Consequently, a modified cross-beam geometry was used (Fig. 

34A ) . This beam arrangement was sufficient to avoid contributions to the 
cross-correlogram from the flow field near both windows. Only disturbances 
in the center of the wake and near the middle of the test section were common 
to both signals. The distance between the beams traversed by the common 
disturbances was 0. 516 inch. 


84 


Flow 


ri2.5v 

\ 


\\ ^ 


Thin Plate 




Figure 34A. Schematic of cross-beam geometry. 
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Figure 34B. Cross-correlogram of local signal in the 
supersonic turbulent wake of the thin-plate model. 




The cross-correlogram computed for this run is shown in Figure 34B. 
The maximum correlation corresponds to a time lag, , of 0. 036 ms 

(0. 034 ms plus 0. 002 ms compensation for parallax). Thus, the most- 
probable transit speed of a disturbance, averaged over the transit distance, 
is 


<U> 



0. 043 
36 X 10"® 


1194 fps 


This speed is 71. 9 percent of the free-stream speed, ; i.e. , 


<U> 

U 


1194 

1660 


0.719 


which is not unreasonable. 

The beams were not crossed in the vertical plane because it was 
desired that the distance between the beams vary across the test section. 
This provided means by which correlations near the windows, if they would 
have been present, could have easily been identified by their location on the 
cross-correlogram. 

This cross-beam measurement was apparently successful from a 
feasibility viewpoint. From these results, it appears that the application of 
the cross-beam method in measuring certain properties of supersonic 
turbulent flow is feasible. 


5. EXPERIMENTAL RESULTS USING A SHADOW-CORRELATION 

METHOD 

In the previous paragraph, it was shown that the interaction zone near 
the test section windows had a significant influence upon the correlograms 
computed from signals retrieved with parallel laser beams. The reasons 
for this are (1) the interaction zone is a zone of more intense turbulence 
than is the turbulence in the plate boundary layer, and (2) the measurements 
were made with the laser schlieren system, which is sensitive to the first 
derivative of the density gradient component normal to the Poynting vectors 
of the beams. 
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There are at least four alternatives offering some relief to this 
problem: 

(a) To redesign the experiment such that the interaction zone is 
eliminated or the intensity of the turbulence in it is reduced to an acceptable 
level. 


(b) To employ a sensing mode which is less sensitive to the inter- 
action zone. 

(c) To use crossed beams instead of parallel beams. 

(d) To study the turbulence in the interaction zone itself, since it 
can already be measured very well. 

Because alternative (a) above is a matter of model or facility design, no 
discussion will be presented here. Alternative (c) has already been proven 
to be a successful solution (Section III, Paragraph D.4). Alternative (d) 
provides interesting material for discussion but falls outside the scope of 
this report. Alternative (c) is of direct importance to the purpose here, but 
will be discussed briefly because of the limitations of the instrumentation 
available for investigation of the shadow-correlation sensing mode. 

In Section III, Paragraph A. 1. , the laser schlieren principle was 
described. It was shown that the ac-coupled time history of a beam was 
directly proportional to its deflection. This was accomplished by placing 
a knife-edge in the beam (Fig. 7). As each disturbance intersects the beam, 
it causes the beam to be deflected in a particular direction, and the knife- 
edge is sensitive to the x-component of this deflection. Thus, it is evident 
that the detailed composition of a disturbance determines the magnitude and 
direction of a deflection. If this is actually the case, then it should be 
possible to sense the changes in composition by monitoring the intensity 
fluctuations inside the beam itself. To do this, the knife-edge in Figure 7 
was replaced with a thin plate having a small pinhole. The diameter of the 
pinhole was smaller than that of the laser beam. The plate is placed 
perpendicular to the beam such that the centerlines of the beam and pinhole 
are collinear (Fig. 35). Because the pinhole diameter is sufficiently small, 
the photodetector is never exposed to the circumference of the laser beam. 
Therefore, the fluctuations in total power monitored by the photodetector are 
a result of the fluctuations in intensity near the circumference of the pinhole 
produced by the disturbances inside the beam. 


87 



Thin Plate 



Figure 35. Schematic of a laser shadow-correlation 
system (side view). 

Two runs were made using the shadow-correlation mode. The beam 
geometry for the first of these is shown in Figure 36A. This run was made 
with the beams crossed in the horizontal plane in the turbulent boundary 
layer on the plate model similar to those described in Section III, Paragraph 
D, 3. The purpose of this run was to compute the contribution to the cross- 
correlogram resulting from signals retrieved from the interaction zones 
near the windows. The cross-correlogram is shown in Figure 36B. At a 
time lag of 0. 160 ms, the contribution of one interaction zone can be detected. 
However, by comparing this correlogram with that of Figure 30A, it can be 
seen that a significant reduction has been achieved. 

The purpose of the second run was to compute the correlogram from 
signals retrieved with parallel beams in the turbulent boundary layer for a 
qualitative comparison with Figure 36B. The beam geometry for this run 
is shown in Figure 37A and the correlogram is shown in Figure 37B. 

Comparing Figures 36B and 37B, it can be concluded that the shadow- 
correlation mode offers a considerable advantage over the schlieren mode 
where the beams must pass through regions of turbulence near or on the 
test section windows. However, the magnitude of the correlation is reduced 
for the shadow-correlogram caused primarily by a decrease in signal-to- 
noise ratio. 

These conclusions must be considered as preliminary because of the 
limitations of the instrumentation and hardware available. 


IV. CONCLUSIONS 


Theoretical discussions and qualitative experimental results have 
been presented in an attempt to assess the feasibility of remotely retrieving 
statistically correlated signals from supersonic aerodynamic turbulence that 


88 






Figure 36A. Schematic of beam geometry. 
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Figure 36B. Shadow-correlogram of positive* time-lag range 
for beam geometry shown in Figure 36A. 
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Figure 37A. Schematic of parallel beam 
geometry (plan view). 
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Figure 37B. Shadow cross-correlogram for parallel beam 
geometry shown in Figure 37A. 
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could successfully be used to compute auto- and cross-correlograms 
exhibiting accurate, reproducible, and readily identifiable flow- related 
peaks of correlation. With respect to this objective, the qualitative results 
presented in this report indicate that optical remote probing of supersonic 
aerodynamic turbulence using statistical correlation is feasible, without the 
use of tracers, in the BWT of MSFC, and, perhaps, in other facilities that 
exhibit similar low facility- induced noise levels. Other conclusions of this 
investigation are as follows: 

(1) Optical remote sensing with parallel or crossed laser beams can 
yield consistent and reproducible statistical results that are related to 
aerodynamic turbulence, the particular relationships requiring theoretical 
analyses and quantitative experimental verification. 

(2) The use of parallel laser beams for retrieving signals from 
two-dimensional turbulence^ increases the power signal-to-noise ratio 
relative to that obtained with crossed beams. However, direct localized 
measurements require the crossed-beam geometry. 

(3) The laser schlieren system presented herein requires isolation 
from the mechanical, and perhaps acoustical, facility- induced noise. For 
this system, the major portion of the facility- induced noise enters the system 
at the light sources. However, it is possible to filter out most mechanical 
vibrations . 

(4) With respect to the data presented in this report, which were 
obtained with a laser schlieren system in combination with a parallel beam 
geometry, the boundary layer on the windows of the test section did not 
significantly contribute to the computed correlograms. However, the 
interaction zone resulting from the interaction of the window and model 
boundary layers apparently dominated the measurements. This result does 
not invalidate measurements taken during this test, since the test objective 
was primarily to ascertain the feasibility of applying electro-optical remote- 
sensing techniques to making qualitative measurements in a supersonic 
turbulent flow. From this standpoint, the primary test objectives were 
achieved. However, before quantitative measurements are made using 
parallel beams, a two-dimensional flow field with minimized edge effects 
should be sought, possibly through improvement in model design and sensing 
mode. 


9. Two-dimensional time-averaged turbulent properties. 
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(5) The shadow-correlation method significantly reduces the con- 
tributions of the interaction zones, and the crossed-beam geometry success- 
fully avoids them. 

(6) The one-shot auto- and cross-correlation methods provide a 
convenient means for measuring the decay history of turbulent structures in 
two- or three-dimensional fluid flows (for detailed comparison of these 
methods see Section III, Paragraph C.4. ). 

(7) The most important advantages of the one-shot autocorrelation 
method are as follows: 

(a) It represents the only method for measuring the decay 
history of turbulence from one sample of data. 

(b) Because it is necessary to monitor only a single composite 
signal, there is a reduction in instrumentation compared to that required for 
a two-channel cross-correlation system. 

(c) Phase matching errors are eliminated. 

(d) Facility operation is reduced by 400 percent because only 
one run of the facility is required to measure the decay history. 

(e) The amount of data handling, storage, etc. , is reduced. 

(f) The overall data-processing time is reduced. 

(g) This method is the only means by which the decay history 
can be observed visually on-line by displaying the one-shot autocorrelograms 
on an oscilloscope as the data are collected and computed from the same set 
of turbulence data. 

(8) The major disadvantage of the one-shot autocorrelation is that 
the signal-to-noise ratio, taken with respect to a particular time history 
contained in the composite, is smaller than the signal-to-noise ratio of the 
particular time history taken separately. 

(9) The one-shot cross-correlation method has the following major 
advantages over the one-shot autocorrelation: There is a higher signal-to- 
noise ratio, and therefore the required integration time for a one-shot 
autocorrelation is 250 to 300 percent longer than that required for a one-shot 
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cross-correlation. (This statement is, of course, restricted to the 
assumptions in Section III, Paragraph C. 4. ) The one-shot method provides 
a means for computing the cross-correlation peak at zero time delay directly 
without the electrically induced time delay necessary with the one-shot 
autocorrelation. 

(10) The computation time required for the one-shot cross- 
correlation is equal to or less than that required to compute the equivalent 
number of individual cross-correlations. 

(11) The method of induced time delay extends the potential applica- 
tion of the one-shot correlation methods and may be used (a) to avoid peak 
interference on a one-shot correlogram, (b) to identify a correlation peak 
partially lost in correlated noise, and (c) to zone the one-shot auto- or 
cross-correlogram, thereby providing a means for identifying a particular 
pair of signals associated with a particular correlation peak on a one-shot 
auto- or cross-correlogram. 

(12) The laser schlieren system described in Section III, Paragraph 
A. 1. is the most reproducible, convenient, and flow-sensitive system 
presently available for studying remote sensing of aerodynamic turbulence 
in wind tunnels. 

(13) The inverted correlation peak shown in Figure 24 and discussed 
in Section III, Paragraph C.3, represents experimental verification that the 
signals retrived from the flow were in fact schlieren signals and thus related 
to the fluctuating density gradient. 

(14) The shape of the correlogram in the region of a peak correlations 
computed from laser schlieren signals is directly related to the position of the 
knife-edges relative to one another and the flow itself. One explanation of this 
is presented in Section III, Paragraph A. 1. where it is theoretically concluded 
that the correlation between signals is proportional to the time-averaged two- 
point product of the fluctuating density gradient component in the direction of 
flow. For parallel beams, this two-point product is integrated along the 
beams from source to detector. 
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APPENDIX A 


DESCRIPTION OF THE BWT FACILITY, 
MODELS, AND INSTRUMENTATION 


This appendix presents a general discussion of important design 
characteristics of MSFC's 7-inch BWT facility, the two-dimensional boundary 
layer and wedge models, instrumentation for remote sensing supersonic 
flows, the general test procedure, and the analog data reduction equipment 
used for this test series. 

The salient features of this equipment and instrumentation and their 
respective roles in data interpretation are discussed below. 


1. MSFC's 7- by 7- Inch BWT 

This test facility is a supersonic blowdown type of wind tunnel 
(Fig. A-1). Dry air at atmospheric pressures (or atmospheric air at 
approximately 14.4 psia with 1. 6-percent variation) is supplied to a 
7- by 7-inch test section and exhausted to essentially vacuum conditions 
(approximately 20 to 29 in. Hg). Figure A-2 shows a schematic represen- 
tation of this facility. Major dimensions of the tunnel are in Figure A-3. 

Two air dryers are connected in tandem to charge the tanks (Fig. 
A-2). The first takes in air at 3000 ft?/min, while the second takes in this 
amount and recirculates an additional 12 000 ft^/min. Two combined tanks 
equipped with a rubber diaphragm liner store up to 60 000 ft?; however, one 
tank is currently not in the circuit. 

Two Fuller duplex and two-stage rotary pumps provide an overall 
vacuum rate of 8920 ft^/min. This results in a 40-minute pump-down time 
from atmosphere (approx. 14.4 psi) to 0. 15 psi for the six interconnected 
cylindrical steel vacuum tanks. The combined capacity for this vacuum 
storage is 42 000 ft?. 

Mach number control is provided by removable nozzles (Fig. A-3) 
machined from solid brass to a tolerance of ± 0. 001 inch. A hydraulically 
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1. Dryer, Type CH 

2. Dryer, Type CHR 

3. Dry Air Storage Tank 

4. Settling Chamber 


5. Upstream Gate Valve 

6. Test Section 

7. Diffuser Section 

8. Downstream Gate Valve 


9. Fixed Diffuser and Duct 

10. Vacuum Tank 

11. Vacuum Pump 




Figure A-2. Schematic drawing of MSFC'S 7- by 7-inch BWT facility. 






PARTICULAR NOZZLE 


A 




activated mechanism clamps and aligns the nozzle in position. The test 
section formed by the nozzle is 7. 029 inches wide and about 7 inches high 
(Fig. A-4). Although several nozzles providing Mach numbers of 1. 54, 

1.99, 2.44, 3,00, 3.26, 4.00, and 5, 00 are available, only Mach number 
1. 99 was used in these tests. The test rhombus, shock angle, mass flow, 
etc. , for this nozzle combination are discussed in Paragraph 2. of this 
appendix. 

The test section side walls which enclose the windows are hinged at 
the top and are activated by (Fig. A-1) hydraulic cylinders. These doors, 
one on either side, extend the full length of the test section. The 6- by 6- 
inch glass windows, 1-1/16 inches thick, are mounted in the doors in such 
a manner that discontinuities between frame and glass surface inside the test 
section (Fig. A-4) are avoided. These windows provide schlieren and 
shadowgraph visualization of the flow region, as well as to facilitate access 
for remote sensing by lasers, ultraviolet, or other light sources. 

Pneumatically inflated seals are used along the nozzle contour 
(Fig. A-4) and at the rim of the doors. Tunnel operation is automatically 
blocked if all seals are not activated. 

Two primary valves (Figs. A-2 and A-3) separate the test section 
from the air supply and vacuum tanks, respectively. Both valves, which 
are hydraulically activated, provide positive sealing when closed. When in 
the open position, a box-shaped extension on the gate is aligned with the 
tunnel duct in such a manner that a smooth continuity of aerodynamic surfaces 
is provided. 

Starting and stopping processes require about 4 seconds. The air 
flow is controlled by the upstream gate valve. This has two advantages: 

(1) starting loads on the model and test section are reduced to a minimum, 
and (2) the tunnel is evacuated before the flow is established, and therefore 
settling time is reduced for the flow since test section pressures are nearer 
vacuum than atmospheric. 

An adjustable outlet diffuser (Fig. A-3) of high aerodynamic 
efficiency increases tunnel running time. The maximum and minimum 
positions of the diffuser are preset for each nozzle (corresponding to a given 
Mach number). At the start condition, the diffuser is open and is auto- 
matically closed immediately after the flow is established to provide maximum 
pressure recovery. 
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The operation of the wind tunnel is automatic. All valves, gates, and 
diffusers are hydraulic, and may be operated through a series of micro- 
switches and relay circuits by one button. By turning the selector switch 
to the manual position, any of the operations can also be controlled individually, 
if desired. 

Figure A -5 shows the details of the control console for the tunnel, 
along with an optional tape recorder. Notice the simplified controls and 
lighted control diagram midway of the center panel. In this visual-control 
schematic, each control element (valve, diffuser, etc. ) is illuminated as it 
is activated either manually or electronically. The close proximity of the 
tunnel control console can be observed in Figure A-1, 

Additional information on the 7-inch BWT can be obtained from 
Reference 6. 


2. Models 

This section describes the design characteristics of the models used 
in this test series. The model most frequently used in this series is a thin 
two-dimensional plate for generation of turbulent boundary layer and wake 
flows. The second model consists of a two-dimensional wedge for simulation 
of base recirculating flows. 


A, THIN PLATE MODEL 

This model is essentially a two-dimensional turbulence-generating 
plate (Fig. A-6) constructed of aluminum. The leading edge has a 9*32' 
taper which extends back for 1. 5 inches. The plate consists of a 1/4-inch- 
thick flat surface for an additional 4. 5 inches. The overall plate, 17 inches 
long, has a trailing edge of 11 inches with a very slow 1*34' taper. This 
model was mounted parallel to the floor on the horizontal centerline of the 
test section (Fig. A-7). It spans the full width of the tunnel for the first 
5-1/4 inches, but has a 5/32-inch cutout or indentation on either side, 
avoiding full contact with the window of the tunnel (Fig. A -4). The model 
showing surface finish, leading edge, and threaded holes (six 8-32) for 
mounting is shown in Figure A -8. 

During these tests, the models were operated with the Mach 1. 99 
nozzle. The flow field resulting from this combination of nozzle and thin- 
plate model is depicted in Figure 3. 
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Figure A-5. 7- by 7-inch BWT control console 
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Figure A-8, Two-dimensional turbulence generating thin plate model. 



The model was mounted as shown in Figures A -4 and A -5. In this 
position, the sonic throat moves slightly downstream (positive x-direction) 
from the nozzle geometric throat, and dual sonic throats (one above and one 
below) are formed on the model about 1 to 1.5 inches downstream of the 
nozzle throat. It can be seen from the shadowgraph in Figure 3 that no 
shock rhombus is formed on by leading edge or reflection from the walls. 

Assuming the throat location described above, one-dimensional Mach 
number, static pressure, density, and temperature variations were calculated 
for the model-nozzle combination. This information is depicted in Figure 
A-9 which also shows the geometric relationship of the model and nozzle 
contour. The freestream flow speed in the test section is approximately 
1660 fps. This information was used for comparison to convection speeds 
measured by the remote sensing system. 


B. TWO-DIMENSIONAL WEDGE MODEL 

The two-dimensional wedge model for simulation of recirculating base 
flows is presented in Figure A-10. The model, about 7 inches wide, 4 inches 
long, and 0. 75-inch thick, spans the full width of the tunnel. Constructed 
of aluminum, the model is relatively smooth and has a surface finish (about 
125). The leading edge has a 9-degree taper. The wedge is attached to the 
tunnel side walls by four 8-32 screws with the axis mounted along the center- 
line, and the tapered leading edge facing upstream. The flow field generated 
by this model is shown in Figure 1. 

3. Instrumentation 

This section describes the instrumentation used for sensing in this 
test series. 

A typical laser-schlieren wiring schematic for remote-sensing 
instrumentation is shown in Figure A -11. Several variations of this arrange- 
ment are possible, depending on the type of flow problem and data desired 
(e.g. , amount of filtering, amplification, general flow conditions, etc.). 

The basic components as seen from the schematic are, in order, a source 
(or laser), photodetector, photodetector power supply, amplifiers, low and 
high cutoff filter (optional for some cases), oscilloscope for observation of 
raw data signal, analog correlator, and finally an oscilloscope and Polaroid 
camera for viewing and recording the final correlogram of the signals. 
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Figure A-11. Typical wiring diagram for parallel beam 
remote-sensing instrumentation. 
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Figure A- 12 shows part of the instrumentation installed at the 7-inch 
BWT. This figure shows the general arrangement of the light source (placed 
for parallel beam measurements), the traversing and elevating stand, the 
analog correlator, and part of the conventional schlieren system from which 
flow studies can be made independent of the remote-sensing system. Figure 
A -13 shows more details of the electronic equipment used for processing 
the raw signals. This equipment is described in more depth in the following 
paragraphs. 


A. PHOTODIODES 

The photodiodes used for signal detection were EG&G (Edgerton, 
Germmeshausen, & Grier, Inc, ) Model S-D 100. This is a fast response 
light detector which operates in the visible and near-infrared spectrum. 
This silicon photodiode, operating over a wide spectral range, has fast 
response time combined with high sensitivity. Spectral response is from 
0. 35 to 1. 13 microns (near ultraviolet to near infrared). The sensitivity 
is 0. 25 Ma/^w at 0. 9jit. Typical values of rise and fall times are 4 x 10“® 
and 15 x 10“® sec. Noise equivalent power (NEP) is in order of 1 x 10“^^ 
watts/^ycps at 1000 cps. A dark current lower than 0. 2 x 10~® amperes can 
be achieved with a bias of 10 volts. 

The temperature operating range varies from - 65° C to +100°C. 

The window is 0. 12 inch in diameter, constructed of Corning 7052 glass. 
Sensitive area is 0. 11 in.^ (0. 073 cm^) . Typical characteristics of this 
photodiode are summarized in Table A-1. 

During testing, these photodiodes are mounted in a convenient 
housing which has a BNC connector. The photodiode housing mounted on 
Uni-Slide elevating mechanisms are shown in Figure A -14. The adhesive 
tape showing on the photodiode holds safety razor blades in position to 
split the beam, thus creating a schlieren effect (see Section III, Paragraph 
A. 1. ). The photodetectors are mounted on an elevating and traversing 
stand similar to the laser source (Fig. A-12). Figure A-15 shows the 
photodetector and mirror arrangement used for the cross-beam measure- 
ment discussed in Section III, Paragraph D. 1. 


B. DC POWER SUPPLY 

The photodiode power supply consists of a 22. 5-volt battery with an 
adjustable potentiometer to regulate the dc current. Normal operation is 
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Figure A-12. Remote-sensing instrumentation installed in MSEC'S 7- by 7-inch BWT. 








Figure A -13. Remote- sensing instrumentation on-line 












TABLE A-1. PHOTODIODE CHARACTERISTICS 



Minimum 

Typical 

Maximum 

Spectral Response, gm 
( 1 0% points ) 

0. 35 


1. 13 

Sensitivity, 0. 9jJ. 
(jua/Mw) 


0.25 


Rise Time (sec@ 90) 


4 X 10~® 


Fall Time (sec @ 90) 


15 X 10-® 


Linearity, n = 1. 0 
(within 5% over 7 decades) 


r" 


Bias (V ) 

1 


150 

Capacitance (Pf@ 90) 


8 


Saturation Photocurrent 
(Peak amps from photo- 
active area, load < 500 ohms 
at 90 bias ) 



0. 12 

Operating Temperature (° C ) 

-65 


+ 100 


less than one-half unit output with a total power up to 1 mW. Figxire A-16 
shows typical noise voltage as a function of frequency. These components are 
housed in a small metal box (2 by 4 by 4 inches) with appropriate BNC 
connectors attached for convenience of operation. 

Figure A-16 shows that, for a 22.5-volt battery, this photodiode has 
less than 5-percent noise from 500 Hz to large values of frequency. Figure 
A-13 shows the photodiode power supply "on-line," installed in its normal 
location. 


C. AMPLIFIERS 


The amplifiers used for these tests were two (one per beam) Redcor 
Model 500, a two-stage type. The operating power of these amplifiei s 
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Noise (Microvolts) 


LEGEND 
0-0 45 Volt BIAS 
x-x 12 Volt BIAS 
Load = 10® OHM 
Bandwidth = 16 H, 


Frequency (kilohertz) 


Figure A-16. Typical noise voltage as a function of 
frequency for the dc power supply-photodiode 
combination used in these tests. 






115 volts ± 10 percent, 60-cycle ac current at 10 watts. The gain for each 
stage ranges from a minimum of 10 to a maximum of 1000 in 10 steps with 
accuracy of ± 0. 02 percent (dc), linearity ± 0, 01 percent, (dc) and stability 
± 0. 01 percent (dc). Bandwidth at 3 db full scale is 1000-kc maximum (dc) 
and 10 cps minimum (dc). Temperature operating range is 0“C to 50°C. 
Table A-2 presents the noise for the Redcor as quoted with 99. 9-percent 
confidence. 


TABLE A-2. TYPICAL AMPLIFIER NOISE 


Frequency 

Noise 

100 kc 

10 ixV rms 

1 kc 

20 juV peak-to-peak 

300 cps 

14 gV peak-to-peak 

30 cps 

4 gV peak-to-peak 


The output is ± 10 volts, ± 10 mA current for both channels, and the settling 
time for a full-scale step input is 20 jusec to 0. 01 percent of final value at 
100-kc bandwidth. Additional features of this amplifier is an adjustable 
bandwidth, high-input impedance, and a bandwidth that is unaffected by gain 
change. It also features a solid-state chopper. 

Figure A -13 shows a pair of Redcor amplifiers on-line; the corre- 
sponding wiring diagram is shown in Figure A-11. 


D. FILTERS 

An Allison variable filter (Model 2BR) was used on each leg of the 
remote sensing circuitry (Figs. A-11 and A-13). This filter features two 
separate networks: a high-pass double-K section filter (0.20 kc) and a low- 
pass double-K section filter (0 to 20 kc). The low-pass frequency was not 
used for this test; rather, the low pass was limited to 30 kc by the Redcor 
low-pass frequency cutoff or to 100 kc when no low-pass filter was used. 
Each filter has 2 controls: an active band switch which changes the cutoff 
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frequency in octave steps and a multiplier dial which tunes the cutoff 
frequency over one octave. The cutoff frequency is that which is attenuated 
approximately 3 db from the minimum insertion loss. The minimum band 
pass without additional insertion loss is one-third octave. The attenuation 
rate of these filters is about 30 db per octave. Each was pretuned to provide 
a minimum phase shift for signals over the frequency range anticipated in 
these tests. Distortion is about 1 percent for a 10-volt input, and at 1 volt, 
it is less than 0. 10 percent. This instrument is recommended for measure- 
ments to 120 db below 1 volt. 


E. LIGHT SOURCE 

The signal source was supplied by two types of lasers: a Spectra- 
Physics (S. P. ) Model 130 Gas laser and a Quantum Physics (Q. P. ) Model 
LS32 laser. These lasers are illustrated in Figure A- 17 in the position for 
parallel beam remote-sensing in the BWT. Figure A-18 shows the quantum 
physics elevated by a Lab-Jack to obtain cross-beam turbulence flow 
measurements in the wake of the thin plate. 

The operational and design characteristics of these lasers are 
examined in more detail in the following: 

• Spectra- Physics Laser. The Model 130 S. P. laser has a helium- 
neon gas-filled plasma tube 27. 5-cm long with an inside diameter of 2. 5 mm. 
The ends of this tube are terminated with optical, schlieren-free, fused 
silica Brewster's angle windows which result in a plane polarized output. 

o 

Wavelength of this output is in the visible red circa 6328 A. 

The standard instrument has a hemispherical resonator with one 
spherical reflector (30 cm radius) and one planar reflector. The resonator 
output power is 0.75 mV cw (minimum) from the spherical end. Resonator 
reflectors are also of optical quality, schlieren-free, fused silica. Reflectors 
are multilayer dielectric coated for at least 99-percent reflectivity (at desired 
wavelength) and are antireflection coated on back surfaces. 

Beam diameter is 1.4 mm at the exit aperture, and diverges less 
than 0. 7 milliradian (145 seconds of arc). Power input requirements are 
provided by a self contained discharge exciter dc power supply. This requires 
a 60 cps 115 volt input of approximately 90 V-A. The laser mode of operation 
normally used for these tests is TEM^ . 
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• Quantum Physics Laser. The Quantum-Physics gas-filled helium- 
neon laser (model LS32) used in these tests has a beam diameter of 1 mm, a 

o 

power level output of 0.75 mW, and a wavelength of 6328 A. In general, this 
laser is similar to the Spectra-Physics Model 130 previously described. 


F. MISCELLANEOUS LABORATORY EQUIPMENT 

Other items of equipment consisted of two conventional Tecktronix 
Model 502A dual-beam oscilloscopes, one for observing and measuring 
peak-to-peak voltage of the raw signal and one for viewing and recording the 
processed or correlated signal from the Princeton analog correlator. The 
signal was recorded on speed 3000 film by a laboratory Polaroid Land 
Camera. By using a manual shutter control, with an f setting of 16 and time 
exposure, the signal was traced on exposed film from a remote switch on the 
Princeton correlator where it had been temporarily stored in a multiple 
capacitor bank. A Digitex digital-type dc voltmeter was used to optimize 
power output by fine adjustment of the photodetector before each run. This 
equipment is shown in Figures A-11, A-12, and A-13. 


4. Test ProceiJure 

Before each series of runs, the equipment was checked to ensure that 
each instrument was adjusted properly. Also, the band pass was selected 
commensurate with the particular test objectives, and filters were phase- 
matched by use of a signal generator and oscilloscope. The next step in the 
procedure was to align the laser source in the desired geometrical position 
and focus the beam onto the photodiode until the photodetector power output 
was maximized. This power output was approximately half of that available, 
since the remaining half of beam intensity was blocked by the knife-edge. 

Prerun checks were made with and without flow to determine the 
extent of extraneous noise that was correlated in the signals. When ground 
loops and extraneous mechanical vibrations were minimized, the desired 
band pass was set, the gains of the amplifier were adjusted to avoid clipping 
the signal, and the correlator gains were adjusted to optimize the corre- 
latable signal to overcome correlator noise, yet at the same time avoid 
signal clipping. Blinking red lights warned when the correlator was over- 
driven. 


Direct current signals were recorded before each run, and the dc 
meter disconnected from the circuit to avoid ground loops. Atmospheric 
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pressure, ambient temperature, and the inches of vacuum were also recorded 
before each run. The tunnel was activated by the automatic start button 
allowing about 5 seconds for the tunnel to settle to steady-state operation. 

The visual display panel indicated the sequence of valves and outlet diffuser 
positioning. Visual observations of these were made during each run. Run- 
time varied from about 40 to 60 seconds for each case. Data were correlated 
on-line and temporarily stored in the analog correlator. Upon completion of 
the run, the correlation function was displayed on the oscilloscope and 
recorded on a Polaroid photograph. 


5. Date Reduction Equipment 

Data taken during this test series were reduced on a Princeton Model 
100 analog correlator. The primary function of this computer is to solve 
the following relationships: 

1 ^ 

V / fl(t) ■ f2(t - T ) dt 

T ^ oo o 

where C is a cross-covariance of the fluctuating portion of random signals 
which contain flow information. Signals f^ and f 2 , respectively, represent 
information from two beams of electromagnetic radiation. If fj and f 2 
are identical, then the result is an autocorrelation where the following result 
is obtained: 


C(t) = lim Jf(t)-f(t-T)dt 

T — 00 o 

This instrument operates as a hybrid computer to solve either of 
the two integrals for 100 n-points of incremental time delay. At . 

The nth point is approximated by 

t-t' 

V (t') ■ V (t' - nt) Adt 

JD 


C (t) = 
n 


RC 


/ 


RC 
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where the following are true: 

• RC = time constant of the averaging circuit. 

• nAt = time coordinates of the computed points. 

• t' represents past history (i.e. , t' > t). 

The three basic operations performed by the computer are ( 1 ) time 
shifting, (2) multiplication, and (3) integration. 

The computed values are stored in a 100-channel analog memory and 
may be recorded during and after computation on stripe harts, x-y recorder, 
or an oscilloscope-camera combination. The latter method was used to 
record data for this report. Computed accuracy should not exceed 1-percent 
deviation from the idealized function for any of the 100 points calculated. 

Other typical specifications of the Princeton correlator are as 

follows: 

• Useful frequency range = dc to about 250 kHz. 

• Averaging time constant (normally 20 seconds) can be varied 
from 0. 1 to 400 seconds. 

• Blinking red lights warn the operator when the dynamic range of 
the input amplifiers is exceeded. 

• Computed values can be stored in the correlator with a decay 
that will not exceed 300 mV in 10 minutes at 25° C. 

• The compilation error is less than 1 percent of true value. 

• Calibration accuracy = ± 2 percent. 

• Gain = 0. 01 to 5. 0 in 1, 2, 5 sequence (with a variation of less 
than 0. 1 percent/hour) . 

• Zero drift = ± 10 mV/hour. 

• Linearity ^ 1 percent 


n 
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• Delay range = 100 fxsec to 10 seconds in 1, 2, 5 sequence. 

• Delay (time base t) accuracy = ± 1 percent at mid-range and 
± 2 percent at five fastest and three slowest ranges. 
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APPENDIX B 


COMMENTS ON THE EXPERIMENTAL NOISE STUDY 
IN MSEC'S 14-INCH TRISONIC WIND TUNNEL 
(TEST NO. TWT 395) 


The experimental investigation presented in this report originated with 
the objective of conducting remote-sensing cross-beam experiments in 
MSFC^s 14- by 14-inch Trisonic Wind Tunnel (TWT) to demonstrate the 
capability of measuring supersonic turbulence parameters using a laser 
source. The TWT cross-beam tests were conducted in the 14-inch tunnel 
special test section that is basically an annular nozzle with provisions for 
attaching an axisymmetric blunt-base body or a nozzle and vertical and 
horizontal windows for cross-beam access. Subsequent analysis of these 
cross-beam data showed essentially negative results because of an unsteady 
flow field in the TWT facility. 

The major problem of the TWT facility was the almost complete 
domination of the data by facility- induced noise. This appendix was written 
to support the contention that excessive facility-induced mechanical and 
acoustic noise dominated the data and thus precluded fruitful cross-beam 
tests in the facility (at least until the remote-sensing system has been 
further developed). The reasons why the 14-inch tunnel was not satisfactory 
for remotely sensed measurements are discussed in the following paragraphs. 

Figure B-1 shows 57 qualitative correlations of various beams and 
accelerometers that were located both inside and outside the turbulent flow 
field of interest. The first three rows across the top depict the autocorrela- 
tion of the individual horizontal and vertical beams and horizontally and 
vertically oriented accelerometers that were mounted on the outer wall of 
the STS. Sixty percent of these show definite periodic trends that are 
detrimental to the extraction of useful turbulent flow information. 

The four remaining rows across the bottom depict cross-correlation 
of the horizontal and vertical beams and cross-correlation of the single 
beam signals with the signal from the individual accelerometers. This latter 
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Figure B-1. 


Cross-beam test in 


MSFC’s 14-inch TWT. 



device gives a qualitative indication of the relative magnitude and frequency 
of noise common between the flow field and structure. As was the auto- 
correlations, these cross-correlations are dominated by an extraneous 
periodic noise. Attempts to filter these data were not successful because of 
an adverse signal-to-noise ratio. 

Figure B-1, Column (1), shows an attempt to obtain an intensity 
traverse. This should have resulted in a correlation function with a minor 
peak at zero time lag since the beams are crossed in a region of low 
turbulence. However, the peak was masked by a periodic signal of 1380 cps. 
Much of this noise can be attributed to structural vibration. This is indicated 
by the cross-correlation of the accelerometer with each individual beam, 
since definite periodic correlations of significant magnitudes apparently 
exist between them. Also, it may be noticed that these beams were going 
through a plastic-sleeve extension of the nozzle outer wall that probably 
contributed to the correlation. Moving beams downstream 2 inches (out of 
the plastic sleeve) apparently did not improve the situation significantly 
[see Column (2)]. The flow conditions and coordinates of Column (3) are the 
same as Columns (1) and (2) except a dominant peak should have been 
obtained at zero time lag, since the correlation volume is located in a dense 
shear layer. In Column (4), the correlation showed a definite mechanical 
or acoustical vibration between beams because one beam was completely 
removed from the flow and directed alongside of the tunnel wall. Further- 
more, Columns (4), (6), (7), and (8) depict various attempts to soft-mount 
the lasers and stands, but this apparently did not eliminate the vibrations. 

Upon analysis of these data, the following conclusions and recom- 
mendations were made: 

• The flow-field, structure, and surrounding environment in the 
14-inch TWT (STS) facility were dominated by an extraneous periodic noise 
that precluded ascertaining useful data at the time of the test. 

• Attempts at filtering these raw data signals and various methods 
of soft-mounting the lasers and detectors did not alleviate the situation so 
that flow-related correlation could be obtained. 

• Remote-sensing measurements in the 14- by 14-inch TWT might 
be possible after (1) more analysis of the 14-inch TWT noise problem, and 
(2) development of techniques to improve signal-to-noise ratio of the 
remote-sensing system. 

• A different facility was needed with a lower facility-induced 
noise level. 
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I, 


These results essentially led to the selection of the 7- by 7-inch BWT 
facility for the immediate future development of the remote- sensing tool. 

This approach appeared to offer the most economical method of expediting 
development of the remote-sensing system for measuring supersonic 
turbulence. 
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APPENDIX C 


UNKNOWN ANOMALIES IN DATA ACOUISITION SYSTEM 


Unknown anomalies should be expected in any research work; the 
research involves identifying and resolving them. The anomalies in the 
system were (1) low-pass filters set too low, thus blocking significant 
amounts of data; (2) use of attenuation after amplification, which, combined 
with low maximum amplifier input voltage, probably resulted in clipping of 
signals; (3) high-pass filter set too low, thus allowing undesirable flow 
noise and facility-induced noise in the signals; (4) filtering low frequencies 
(e.g. , 60 cps) with large amplitudes before the run, probably causing 
clipping of signals during the run, which was not detected after the clipped 
signal was filtered; and (5) photodetectors and other possible points of 
contact between grounds were not isolated (not always a source of trouble). 

As can be seen, the problem areas were associated with the ampli- 
fiers (i.e., clipping) and filtering methods. 
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